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GHAPIER - 1 


INmODUaCIQU 

Ibe present tliesis makes contri'buticns to the following 
five areas of mathematical hiosciences • 

(a) Mathematical Models for Population Growth 

(b) Difference Eq.uatians Models in Ecology and Epidemics 

(c) Relative Stability of Difference and Differential 
Equation Systems 

(d) Application of Bifurcation Iheory to Population Models 

(e) transition to Stochastic Models from Deterministic 
Models . 

these contributions are given in chapter IL-VI of the 
present thesis. 

In Chapter II, we give more than half a dozen models 
for growth of populations of microorganisms. A microbe or 
a microorganism is too small to be seen by the naked eye and 
is less than 0.1 mm of diameter. However inspite of their 
small size, these microbes play an iii|)ortant role in the 
following disciplines: 

( 1) Per mentation technology where microorganisms are 
used for producing foods, beverages, antibiotics, vitamins, 
plant growth regulators, flavour enhancing compounds, amino 
acids, enz3n3ies, polysaccharides, single cell proteins and so on. 
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(2 ) Sanitary and Envjjccnmental Engineer lag; wliere 
micr ocjrganisms ane used to reduce v/ater p'ollution, 

(3) Ore and ITuel Prooessiag where microorganisms can be 
used to leach certain, undesirable substances from ores e4g, 
they can be used to remove sulphur from coal. 

(4) Bioconversion of Solar Energy where microorganisms 
can be used to first absorb solar energy and then they can 
be used as a fuel. 


Our models are based on the esperimen tally determined nature 
of the function k(c) where lc(c) is the specific growth 
function for microorganism or i.e. it is growth rate per 
microorganism and c is the cancentraticai of the substrate 
from which the microbe derives the nourishment. 

In the literature various forms are given for k(c). 

(a) Monod’s function 1949 ^ 

(b) leisser's Bunctioa [1 1942]] 

(c) Jost et al Function []]19’73']] 

(d) Shehata and Mar r Function ^1971]] 

(e) Babes et al Function C^^VS]] 

(f) Moser's Function [] 1958]] . 

Corresponding to each of these, we get a mathematical model 
to determine the population size ]Sf(t) as a functicai of 
time, Yfe find that in each case, the population has a 
limiting size and each grovrth curve has a point of inflexion 



Since most of the growth curves observed on practice have 
these features, the large variety of models gives us a good 
scope for fitting of data for different microbes and different 
substrates, ¥e have also investigated general theoretical 
models which give rise to sigmoid growth curves with points 
of inflexion and limiting population sizes. 

All the models in the second chapter are in terms of 
differential equations. However difference equation models 
are also used in population dynamics, ecology, epidemics etC', 
These are relatively sirpler to understand and manupulate 
because these do not require a knowledge of calculus. These 
models can be under sta&d even by non-mat hematic ally oriented 
demographers, ecologists and ep idimologists and can enable a3J. 
college students to understand the great problems in this 
discipline. Both these types of models have been used in 
the literature and both can represent reality. In fact the 
validity of a model can only be tested by conparing its 
predictions with actual observations. Moreover differential 
equations are usually obtained by first writing difference 
equaticms and then using a limiting process, Bor numerical 
integration we have again to transform the differential 
equations to difference equations. In recent years there has 
been an increased interest in modelling ecosystem in terms 
of difference equations n^^Ihis interest stems from the 
recognition that many species have in^ortant discrete components 



in their life histories, such as synchronized hreedang, 
seasons and seasonal death. Difference eq^uation models 
are more or less analc^ous to the more common differential 
eq- nation models. 

In Chapter III we accordingly investigate the following 
difference equations models. 

1. Growth of a population 

2, Growth of a populaticxn with limited resources 
3'. Influence of pollution on population growth 

4. Influence of age structure 

5. An alternative method of taking age structure Into account 

6. Prey-Predator and Host-Parasite models 

7. Interaction between predators and prey 

8. Leslie’s equations 

9. Gcmpetition between two species 
10« Oonpetition among three species 

11. Epidemic model v/ithout removal 

12. Epidemic model with removal 

13. Growth of population with time-.delay 

14. Growth of population with time-delay with limited resource 

15. Prey-Predator model with time-delay. 

In each of these we have drawn graphs and discussed important 
conclusions. 

One of the important problems in ecology is the stability 
of the system. While both difference and differential 



equations can be used to model ecosystems > these models may 
not produce identical stability behaviour for corresponding 
models. Corresponding to difference equation model, we can 
try to obtain a differential equation model, by using a 
limiting process and corresponding to each differential 
equation model we may try to obtain a difference equaticaa 
model by converting differential equation to difference 
equation. In Chapter IT, we discuss the relative stability 
of some correspondence models, For defining the analogous 
discrete system, May !>■??] has replaced the derivative by a 
forward difference operator with step size h > 0 so that 


dN.(t) 
.. . 3 . ..,1 
dt 


'Vt ~ 

h 


(T-1) H. 


where 

!M^(t) = ]J^(t+h) 

and time has been normalised so that h = 1 May 0-973 has 
stated that: 


'’It is widely understood that difference equations 
tend to be less stable than their differential 
equation twins, because the finite time lapse between 
generations of growth will have the destabilizing 
effects associated with any time in an interacting 
system. Our discussicm makes it explicit, clearly 
stability of the difference equations system implies 
stability of the differential equation one, but the 
converse is not necessarily true" • 
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However Driessche , , using the cperatianal rule 

■ h 

• (I-l)N. (1+1) 

has shown that the resulting analogous discrete system is as 
stable as the differential equation system. 

In this chapter we examine the operational rule 


dtl.. 


4 pfi ST 


where (p ^ -1) is a parameter, to compare the stability of 
the continuous and discrete models. Ihe results of 
Driessche [l97^ and May Pl 97^ should follow as particular 
oases of ours and we may expect to get additional insight by 
using other values of p , We have shown that the relative 
stability of the difference and differential equation models 
depends on the operator used and we have demonstrated a class 
of operators for v/hioh difference equaticai model is more 
stable than the differential equation twin, This discussion 
illustrated with an age-structured population model with 
three age-groips, Ihese three age groups are the pre- 
reproductive, the reproductive and tte. post-reproductive groups. 

In Chapter T bifurcation theory is used to discuss the 
stability of some biological models. Bifurcation theory made 
significant progress and has found applications in various 
established fields and to many new engineering sciences. 
Bifurcation theory has been widely acclaimed and has found 
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fruitful application in theoretical biology* It is of . 
particular interest in bifurcation theory to study how the 
solution and their multiplicity changes as a parameter varies, 
nature of the equilibrium point changes when parameter values 
changes slightly in this case* For two dimensional case if 
the roots are V+io) , we get an unstable focus if V > Oj 
a centre if V = o and a stable node if Y < o. V is a 
function of the parameters and for some values of the parameter, 
Y can be zero, A slight change in the value of the parameters 
would change the nature of the equilibrium. In this chepter, 
a set of ccnditions for a pair of eigenvalues of a stability 
problem matrix to be purely imaginary are obtained. Ihese 
conditians are applied to discuss the existence of bifurcaticn 
point for (i) time-delay prey-predator model and (ii) growth 
of a single population under the effect of pollution. 

Just as there is a correspondence between differential 
equation and difference equation modelA^there is also corres- 
pondence between deterministic and stochastic models. 
Corresponding to a given deterministic model in terms of 
differential equation, we can find a stochastic model. In 
general we can get a system of differential-difference equation 
for this stochastic model and can solve the problems by first 
solving a partial differential equation for the generating 
function. liYe find a number of such partial differential 
equations in the siocth chapter but unfortunately most of these 
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appear 1;o be intractable either because these eq:Uations are 
non-linear or because the initial and boundary conditions 
ar« not sufficient cnr because the p-e3?tial differential 
equation involve unknoTOi probabilities. Here the partial 
differential equation have been used to determine some 
interesting relations between the e^qjected values of the 
variables in some special cases for the some of the models 
discussed, 

She work done in the present thesis has been published 
or accepted for publication or communicated for publication in, 
the form of the following papers* 

( 1) Some Mathematical Models Por Population Growth 

published in Indian J. Pure Appl. Math., Yol. 10 » 
Ho. 3, p, 277-286, March 1979. 

(2) Bifurcation Theory I'or Tvro Population Models 

published m Indian J. Pure Appl. Math., Yol. 9, 
lo. 8 (1978), pp, 787-796, 

(S) Relative Stability of Difference And Differential 
Equations System published in Rational Academy 
of Science Yol. 2, Ho, 4, April 1979. 

(4) Difference Equation Models in Ecology and i^idemics 

(To be communicated), 

(5) Some Mathematical Models of Population Growth ~ II 

(To be communicated). 

(6) Stability of Differential Equation and Difference 

Equation Syston (Communicated for publication). 



CHAPTER - 2 


MATHEaJIATIGilL MOIELS FOR POPULATION GROV/TH 

2 *1 INTRODUCTION 

Among tlB many aspects of a populaticn, ecologists are 
concerned with its grovrth taking into account its fertility^ 
mortality, gra^rth form and regulation (fluctuations, oscillationsj 
dispersal etc.). They are also interested in inter and 
interaspecific interactions among populations ( cooperation, 
predation, conpetition etc,). Since its inception, this 
field of population dynamics has had a strong mathematical 
orientation in formulating and testing descriptive and 
predictive models. 

It was observed by Malthus [jL’79^ many years ago that 
the human population tends to grow in geometrical progression 
or the rate of population increase is proportional to the 
population itself. 

By 1830' s it was noted by several persons that the strength 
of ejsponentiat ion process seems to-be declining. Various 
hypotheses concerning the resistance or the obstacles put in 
tl^ way of the indefinite development of the population can be 
made. Quetlet suggested that saturation might even be expected, 
Verhulst [l83^ produced a correction factor to the Malthus 
eq^uation which v/ould inply saturation 
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^ = KH [>(u/e ) n , (1) 

where 0 is the saturation population* Of course as 
0 ^ 00 this reduces to Malthus equation. 

The first one to atteicpt an adequate mathematical 
representation of the normal rate of growth of the population 
of united states was Pritchett. Taking the census data from 
1790 to 1880, Pritchett fitted by the method of least squares 
the following equation 

P = A + Bt + Gt^ + Bt^ (2) 

where P represents the population and t the time from 
some assumed epoch. Pritchett got a very accurate represen- 
tation of the population between the dates covered. As 
will presently appear, this curve did not give, even v/ithin 
the period covered as accurate results as a really satisfactory 
curve should have done • 

During the early 1920’ s Pearl and Reed Q-924] analyzed 
the population growth of many European countries by en^loying 
the form 

net) = 0^+ C^/ [1 + Og (2) 

The extra parameter was introduced so that the two 

growth regims could be characterized by one equation. G^ 
is the value of N(t) as t 
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A family of emperical curves which saturates either more 
slowly or more rapidly than the initial form of the rise can 
he derived from the followring general izaticxn of the Yerhulst 
e quat ion DS3C 

g=KH Cl- (4) 

which is the same as (1) when a = 1, As a *♦ o> (4) gives 

^=-Olog (N/e) (5) 

which is known in the populaticaa literature as the (^cmpertz 
p.82^ equation. This equation never becomes popular 
because this was more difficult to fit by least squares to 
census data than was the Yerhulst equation* 

In view of fluctuating population driving forces, it 
would seem that a more appropriate model for population 
variation would be characterized by the equation 

^ = HT G(N/e) + N P(t) (6) 

where F(t) is considered to be a random function of time , 
which reflects changes in the economy and altitudes. 

One of the most successful models for e 3 q)laining the 
growth of populations of bacteria and even of humans is the' 
so-called logistic model (1)» Its characteristic features are: 

(i) a limiting population size and ultimate zero growth rate 

(ii) an S-shaped graph with a point of inflexion. 
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Many observed populations show these features but do not 
fit the logistic model very well. Smith p.96^ gave 
another model with these features giving a better fit in 
some situations. In the present chapter half a dozen new 
models for microbial population growth are given, Ihey 
are based on the esperimen tally determined laws for the 
dependence of the specific growth rate of a population on 
the concentration of rate-limiting substrate . It is shewn 
that though these laws are quite different, they all lead to 
logistic- type growth curves with points of inflexion and 
1 imit ing p op ulat ion si ze s , 

2.2 BgORIMGS OF MOBBLS OF PaPUMTIOH GROWIH 

Actuaries and demographers are interested in models of 
growth for human population for predicting e 2 {pected duration 
of life at various ages and for estimating future population 
trends (Pollard jjL97^ ) . Bioeconomists are interested in 
models of growth for populations of sheep, fish, forests and 
other renewable resources for the saJee of their optimal 
esploitation (Olark p.97S] ) « Medical scientists and biologists 
are interested in models of growth of bacterial populations 
for the sake of controlling diseases, and epidemics and for 
genetical studies, Chemical engineers are interested in 
models of growth of populations of all microorganisms for the 
role they play or can play in four major areas viz. 
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(i) [Fermentation technology where microorganisms are used 
for producing foods, beverages, antibiotics, vitamins, plant 
growth regulators, flavour enhancing conpounds, amino acids, 
enzymes, polysaccharides, single-cell proteins and so on. 

(ii) Sanitary and environmental engineering where microorganisms 
are used to reduce water pollution. 

(iii) Ore and fuel processing where microorganisms can be used 
to leach certain undesirable substances from ores e.g. they 
can be used to remove sulphur from coal 

(iv) Bioconversion of solar energy where microorganisms can 
be used to firsb absorb it and they can be used as fuels . 
(Erederickson and Isuchiya [197‘3 ) • 

One specific example is the release of methane gas from 
bio-gas plant which is due to the action of microorganisms. 

Eor increasing the amount of gas released, it is necessary 
to understand the mutual relationships between the growth of 
microorganisms, the change in the substrate and the amount of 
gas produced. 


2 ,3 THB LOGISIIO OURVE 

The simplest growth model is 


ws = 


aE[, 



where l!l(t) is the population at time t and a is the excess 
of births over deaths per unit time per individual. 


This 
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is known as law of Malthas and gives an integration 

ll(t) = (8) 

This model is valid when there are unlimited resources 
available to the populaticn. However limitation of 
resources suggest the intr oducticai of an inhibiting term 
on the right hand side of (7) to take into account the 
competition for resources among members of the population* 
The sinplest modified law is the logistic law 

M = aH - bU^ (9) 


Integration of (9) gives 


^e ^ r ^-at 




l+(^ *-l) ®* 


a 

ts if*-* 

e D 


( 10 ) 


The logistic curve is S-shaped and Ims a point of inflexion 
¥/hen 




( 11 ) 


curve is shown in figure ^2,1 . 


This law explains wery well the growth of a bacterial 
colony in a nutrient medium or bifetfa, provided the initial 
nooulum is taken from an exponentially growing colony 

Gauss &93|1 . Rubin ow [j.975] • It has also been successfully 
used to fit data for human pqjulaticais. Equation (9) is 
sometimes replaced by 

^ = a(t)H - b(t) 


(IE) 
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v/bioli gives on integraticEti 

i = ^ b{-fc) at3 (13) 

where ^ 

a(t) = / a(t)dt (14) 

0 


2 .4 MOBSL OF StlllH 

According to (9), ^ ^ should he a linear function 
of IT . Smith [1^,63] and Pielnu [l96^ found in his experiments 
on bacterial cultures that while the graph of N against t 
was a sigmoid or S~shaped curve, the graph of ^ ^ against 

1 was not a straight line hut it was a concave curve • Ife 
argued that the term a-hN in (9) should he replaced hy the 
rate of food supply not being used by the populaticm i.e. 
it should he replaced hy y(S-P)/S where P is the rate at 
which food is being used and S is the saturation rate. He 
assumed that 


P = XN + 


V 


dt 


so that 


1 W ^S-P 


Y( XK - xn ~ ) 

XK 


(15) 


(16) 


or 


1 m Y(K-ir) 


(17) 


where 
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\> ts , This model gave a good fit for his data. 

If Nq < K, we find that ^ ^ Oj the population goes on 
increashig at a decreasing rate and ultimately as K -* K, 
the growth rate -► 0 so that K is the limiting size of the 
population. If > K, the population decreases to the 

limiting size K, 


Also 


1 

Y 


2 

d^N 


k2 


(18) 


_ ^ V 

dt^ (K+v^ 
so that there is a point of inflexion when 


1= CT “ H 

As 0» ^ ^ s-s V -* ooy ^ 0, so that the point 

of inflexion always oocurs wlBn N lies between 0 and K/2. 

As v.(or u/x ) increases, the point of inflexion occurs 
for lower and lower values of IT and if 

V > ^ ^ (20) 

there is no point of inflexion. 

Thus Smith’s model is lifcely to fit the data better v/hen 
the point of inflexion on the sigmoid curve occurs before half 
tlB final population size is reached. We may also note that 
(17) has one parameter more than (9) and in fact includes (9) 
as a special case . 
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2.5 MOIEIS BASED OS MONOD^s HIPOIHBSIS 
The law of growth is given by 

^ = K((3)H (El) 

where K(0) is the specific growth rate and 0 is the 
concentration of the nutrient in the medium. Monod 
pL949>5(^and E-ubinow [19V^ assumed that the mathematical 
form for K(G) was expressed by the Michaelis-Menten 
equation of en 2 yme kinetics i.e. 

K(0) = (EE) 

The constant is the maximum specific growth rate of 
the population, obtained v/hen 0 » K,, and the constant K 
is so-called Michaelis constant. K is the concentration 
of rate-limiting substrate at which K(G) is one half of 

Kjjj. Many authors have found that Monod’ s law (22) does not 
fit growth rate data very well and they advanced other laws 
that fit data better. 

2 .6 MOIiBl BASED OS THE LAW OP TEISSEB. 

Teisser [1941] found that equation 

g(g) = \ C (- ° ^ ^ ) J (23) 

fitted his data quite well. 

The law of growth in this case is given by 




( 84 ) 
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She caticentraticai G Is reduced due to the consunptiou 
of the nutrient by the bacteria according to the law 


|2 = _ i K(C) H (26) 

SO that consurption rate is being assumed proportional to 
both the population size of the bacteria and saturation 
function of the nutrient. Ifere y, the yield coefficient 
is the amount of biomass produced per unit amount of rate- 
limiting substrate consumed* 

Prom equation (21) and (25) 



which gives cm integration 
H ta y(G^-G) 

where is the constant of integration. 

Prom equaticm (24) and (27) we get 


dt 


= K^Oexp {-(Cj^-H/y) } ;;3 B 


(26) 

(27) 


(28) 


After integrating equation (28), we find 
X dX 

/ ■ “ - T (29) 

^oXCl-«sp {~B(1-X)}2] 

where 

G.y « 0, X = Xy, B = , t == jr t, X= § (30) 

i X 0 

and X^ is the initial populaticm. 
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Pig, 2,2 shov/s the population curves for different 
values of B when the initial population = 0.1. 

A point of inflexion, which is the characteristic 
feature for the population growth can be found out by 
d% 

epLuating equal to zero, 

dt^ 

Differentiating once the equation (21) we get 

G-fK.(0)I^K(0)3 (31) 

Por point of inflexion 

yK(0)=mC‘(G) (32) 

with the help of equation (23) and (32) we get the point 
of inflexion X* will exist if 


1-(BX*+1) = 0 

Applying condition (33), points of inflexion in the sub 
figures of figure 2.2 are given below 


Sub figures 

B 

Point of inflexion X'^ 

1 

0.30 

.518 

2 

0.25 

.515 

3 

0 . 2 o 

.512 

4 

0.15 

,564 

5 

0.10 

.506 
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2 

Yfben B is very small (33) gives, on neglecting B and higher 
pov/ers 

1-2X^-B jj;*^ + I “ 0 

so that as B”*0, X ^ 

i*e* for very small values of B, point of inflexicxn will 
occur near half of the final populaticai size* 

Let 

f(X) == 1-(1+BX) (35) 

v/e give different values of X and find the sign of this 
function 

f(0) = l-e“^ > 0 


f(l) = -B < 0 


nh 


L.(l+ iB)e 


B 


B 


?B 


= e 


[e^ -(1+ ± B)3> 0 


(36) 


These show that the point of inflexion X* lies between 
1 

and g> 
so that 

^ < X^^ 1 

As B -♦ «>, X ”♦ 1 


g}_h us a point of inflex i on always exists and occur when t he 
population size is more than half the final population size . 



23 


2,7 MOjEL based on THE DAW OE JOST ET Al 


Jost et al |JL97^ obtained from their experiments 
the law 2 

^ (Kp+C) 


(37) 


Growth of the population in this case is given by 

ffl ^ 

(Kj+'iO'CKgVcrj ® 


(38) 


Erom equation (27) and (38), after integration we get, 


/ IfcgJlgrM ax = T (39) 

Xq C^X(1-X)2 

where „ 

O^y = 0 

H = E^y + G, G = K^y + 0, X = r ,T=K^'fc (40) 

c 

and Xq is the initial population 

Eigure 2,3 shows the population growth for different values 
of H and the initial population X^ = *02, G = 3,0, 0 = 1,0. 

Values of H in sub figures (1), (2), (3), (4) and (5) 
are 2.4, 2,6, 2,8, 3.0 and 3.2 respectively. 

Erom equation (32) and (37) v/e find that point of inflexion 
will exist if 

•tCX) = (IL.®)(G-0X)(1-3X)+^ (G+H;-20X)(i-X)=0 (41) 

To find the point of inflexion, we give different values of X 
and see where the sign of the function change 
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IPor X = I 

n|)> 0 

and for X = 1 

( 1 ) < 0 

As such a point of ±nflexiaa always exist and ocours, after 
o ne-. third the final population size is reached 

If G = H, (41) gives 

(1) = (G-OX) [2 (G-OX) ( 1-3X)+2GX( 1-X) ]= 0 (42 ) 

Since G > 0, the first factor camot he zero when 0 < S < 1 
end therefore the point of inflexicxn is given hy 

(G-X)(a^3X) + 2XX1-S) = 0 

or 

X^ + X(1-.3G) + G = 0 (45) 

G* 

where G « •“ 

0 

Equation (43) has two positive roots, one greater' than and 
the other less than unity, She latter root gives the point 
of inflexion. So that 

X'^ * I C (3G-.1) - f (9G-.1)(G-1)3' (44) 

For some different values of G, the points of inflexion are 
G = - i I z 

X* 1 I .44 

and as G ^ X* ^ • 
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Equation (44) determines the point of inflexion for the law 

Eor the general lav/ (37), the hasio equation (41) is a cubic 

in X and its root lying between ig- and 1 can be determined 

o 

numerically. Using equation (44), the points of inflexion 
in sub figures of figure ,2.3 are given below 


Sub figures 

H 

Point of inflexion X* 

1 

2.4 

.4045 

2 

2.6 

#4005 

3 

2,8 

.397214 

4 

3.0 

.3944 

5 

3,2 

♦3920 

MOIEL BASED ON THE 

lATf Of SHBHATA AND mBR 


Shehata and Marr p.97l] found from their experiments 
the law 

or a generalisation of this involving the addition of further 
terms of the same form.. They pointed out that this equation 
would arise if growth involved different but parallel process 
of substrate i^taEe • 


The growth of the population in this case is given by 


dH 

dt 


K^G 

- (.A— 

'“K^+C 


+ 



(4V) 
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From eq^uaticin. (27) and (4:7) after integration we get 

/ (P-gX)(E-01) ^ ^ 

X(l-X)[B-A:g 

where D =s K.y + 0, E w K^y +G, Xs|., T~ict 

' Q m 

and A and B are positive constants 
Figure 2.4 shcjws curves for 


.02, 0 = 1.0, D = 2,2, E =s 2,4, Am 5,0 and 
B M 18, 16, 14, 12, 10 


‘fhe points of inflexion are given hy 

(X) = K^(iLx)^(2xE-5-.X^) + X2(5 uX)^( 2EX-E-X^) a 0 (49) 

where 5=2, E = !• 

0 0 


So find the point of inflexion, we give different values to X 
and see when sign of the function change 


for 



and 


1 


V'(O) « -K^E^ 5 - E < 0 
,(|)=.'^( E -|)".^( D .|)^<0 

(50) 


4'(1) a K^(E~ 1)^(5-.!) + K2(D-1)^(E-1) > 0 
since 5 > 1, B > 1 

As such a point of inflexion always exists and occur when t he 
population has reached mo r e then half the population size . 

In fact ''l»(X) is negative v^hen X<^, 'l'(l)>0, t|'(“) < 0 

and as such ’t'(X) = 0 has cne positive root hetvfeen ^ and 1 
and at least cne positive root Between 1 and *», 
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If D = E, then the root of (49) are 


X = 5, 5 > 1, 


X = D + 


Y d'^-5 


(51) 


The point of inflexion occurs at 

X* = 5 (52) 

hut in this case the lawxShehata and Marr reduces to 
Monod’s law 

Eor the Eigure 2.4, the points of inflexion are gi-ven hy 


Sub Pigures 

B 

Points of inflexion X* 

1 

10 .0 

,552 

2 

12,0 

.535 

3 

14,0 

,513 

4 

16 ,0 

.507 

5 

18,0 

.505 

MODEL BASED OE 

THE LAW OE 

BABES El AL 


Powell [ 196 pointed out that a diffiusional resistance 
to substrate transfer through the medium around would affect 
the growth rate. When ilond’s model Y/as modified to account 
for this, Powell obtained 


K(G) 


KJK+Ih-G) 

— sir~ 


D- 


fu 


ziissr 

(K+Ih-gT^ 




(53) 


where L is the parameter having the dimension of concentraticai 
which depends on the diffusicnal resistance around the cell, 
Powell’s equation reduces to that of Monod of 1 is set equal 
to zero. 










H 


MMwiiMi 







50 


Recently, Dabes et al p.97^ ha’ve proposed a number of 
nev7 growth rate eq[uaticins 

K(0) = (54) . 

Growth of the population in Powell and Dabes models are 
similar, ifere Labe model is ta]d®n* Prom equation (21), 
(2V) and (54), after integration we get 


/ 'Ti' ' I I II I « f (55) 

xc (H -CX)^V (h',.I 3S)^-G (l^X) 

where 0 = O^y, X « I- , t = K^t, (56) 

G 

(Bfi^)y + 0 * H*, 4,M^Q « G* and is the 

initial population. 


Pigure 2,5 shows the population growth for different values 
of H‘ when the Initial population X^ = .02, G* a 1.2, 0 = 1,0. 

Values of H* for sub figures (1), (2), (3), (4) and (5) 
of Pigure 2,5 are 2*4, 2*6, 2*8, 3*0 and 3*2 respeoti’vely » 
Prom equation (32) and (54) we find that point of 
Inflexion will exist if 


2(H‘-0Xr.l) T(H'-CX)"'^rG*(l-X) + X [50(H*-GX) - g’ ]] 

- 2 [;(H’..CX)2 -,G'(a^X);] « 0 (57) 

using equation (57) we can find the point of inflexion of 
sub figures contained in Pigure 2 ,5 • These are given 
below 








mm 




ip.. b‘0i|6}'ficiCcj 











Sub figures 

H* Point 

of inflexion 

1 

2.4 

,56o4: 

2 

2,6 

.5548 

3 

2.8 

,5502 

4 

3,0 

.5463 

5 

3.2 

,5430 

MODEL BASED ON 

THE LAW OP MOSM 



We give here two more models of the same type based on 
oonsideration of enzyme kinetics and compare these with the 
existing models. 

A more general form of (22) is given by Hill’s equation 
or Moser law Q.95^ 


K(0) = K 

(58) 

Prom (21), (27), (58), we get 


"* K^+(o^-.N/^)^ ■l;^+(a-N)^ 

(59) 

where I = , 0 = O^y 

(60) 

we note that 



i) If < 0, then ^ > 0, N goes on increasing and as 
N "* 0, ^ 0 so that 0 is the limiting size of the 

population. In our model oannot be negative 
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ii) jDhe graph ^ ^ against IT is not a straight line 
hut may he a concave or convex curve • 

iii) There are two eqLUilibrium points viz* N = 0 and 
IT at G The first is unstable and the second is stable. 
If IT is Initially zero, it always remain zero, but if 

is slightly greater than zero, then H grows with time and 
approaches 0 asyn^) totically , 


iv) 


A 

dt^ 


HI 


I (^H)^'^^+K'(^lT)-nKf^ ^ 


(61) 


let 


f(IT) = (0-^)^***^+ K (0 - IT) - nlTK 






^^)n+l ^ (i„x)^^^+A (I-OT: X) } « (0) 


_ n+1 


where 

„ H , _ k'^ k'*' 

5 ’ ^ ^ 

In general X a “ « 1 and we oan aasuiae 


(62) 

♦ (X) 

(65) 

(64) 


*(X.) = (]-0C„)“+iA (l-nJT X„) > nAX, 


(65) 


Also i|>(l) is negative and so the population curve has a 
point of inflexion at X* between X^ and 1. The pooition 
of the point of inflexion will depend on A and n» 


I'or a given n, we can choose A so that the point of 
inflexion is at a given point X'* where < X^ < 1 

v) For a given observed population curve with a point of 



inflex icai and a limiting population size, 0 is given “by tlB 

liniting populaticoi size, then = N^/G is given by the 

initial population size, fbe existence of the point of 

1 

inflexion shows that n > — r - 1 and then A is given by 

X'^ 


A= 

we can choose n to get best fit to the data* 


( 66 ) 


vi) The model is more general than the logistic model as 
well as Smith’s model in the sense that if can give population 
curves with points of inflexion both before and after 0/2 
while logistic model gives a point of inflexion only at "B/E 
and Smiths model gives only before G/2 » This is not 
surprising since our model is a four parameter model, whereas 
the earlier model were two and three parameter models 
respectively, 

vii) Prom (59) and (64) 


1 M = (^x)^ 

^ A+(l^)^ 

where t =s K t 


Integrating 


X 


/ 

X. 


A 


^X(l-X)^ ^ 


if n - 


when n = 1, this gives 


(A+1) log ~ + A log 


(67) 

( 68 ) 

(69) 


=: T 


(70) 
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If n Is a positive integer greater than unity, we get 
(A+1) log A log + ••• 


(n-l)(l-X)' 


n- 


tH; 


X 


= T 


(VI) 


When n is not an integer, we have to integrate numerically. 

We note that t increases with A,X and, n and decreases with X^* 


Figures 2, 6-2, 9 show population growth curves for 
different values of A and n for X^s 0.1, 

From eq^uation (63) we can see that point of inflexion 
will exist if 


<^(X) - (1-.X)^'^^+A (i^n+1 X) = 0 


(V2) 


Points of inflexion of Figures 2 •6-2 ,9 are given below 


A=0 . 

FIG, 

5,Z^=0.1 

[?.g 

A=2. 

FIG 

0,X =0.1 

[5.vl 

A=3.0,X^=0.1 
FIG [2,^ 

A=4 

FIG 

.0,K =0.1 

[2 .a 

n 

Point of 
Inflexicn 

n 

Point of 
Inflexion 

n 

Point of 
Inflexion 

n 

Point of 
Inf lexica 

.33 

•860936 

.33 

,796977 

.33 

,78443 

CO 

CO 

• 

,77737 

cn 

o 

.792602 

,50 

.716881 

,50 

,70269 

',50 

,69477 

,66 

.735166 

.66 

.654111 

.66 

,63935 

.66 

.631168 

1.0 

,633975 

1.0 

,55051 

1.0 

,53589 

1.0 

.527864 

2,0 

.446426 

2.0 

.374183 

2,0 

,36216. 

2,0 

.355629 

3 #0 

,343102 

3 .0 

.283030 

3 ,0 

,27324 

3.0 

.26794 

4,0 

,278310 

4,0 

,227509 

4,0 

.21933 

4,0 

.21491 

5.0 

.234002 

5.0 

.190172 

5.0 

.18316 

5.0 

.17939 
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2,11 MOTHER pgPULAg l On MOIiBL 


Another equation for enzyme moderate reaction is 

K(0) =K 

^ (1+0)^ + BO^- 

Brom (21), (27), (6o), (64) and (73) 

^ ( L-X) XC ( 2H- 1>.X) ^+B( 1-X) ^ ] 

( 2H- l^X)^ + B( 1-X) ^ ■ 


(73) 


(74) 


where 


z = ^ 
0 


(76) 


Equation (74) can he easily integrated by quadrature 
Por the point of inflexion we find 


d^X _ 


’I'(X) ^ 


T=r^ 


(z+l-X)^+B(l-X)^3 


(76) 


where 


iJ)(X) = [; (z+1-X)\b( 1-X)'^3 X f {(zH-l-X)^(l-X) + B(l^X)^} 

- X { (2H-1-X)‘\3(1-X)(2H-1^I)^+ 4B(1-X)^2 

+ X Q(l-X)(&fl-X)^+ B(l^X)^>x{4(&i-l-X)^+4B(l-X)^>I] 

(77) 

so that 

^(1) = - z*^ < 0 (78) 


Also ’l'(X|^) >0 if X^ is sufficiently small. As such 
the population curve is an S-shaped curve with a limiting 
population size and each curve has in general a point of 
inflexion • 



¥e have got here a t hree—parameter family of curves 

Figures 2.lo~2.13 show populations curves for different 
values of B and Z. Points of inflexion in these figures will 
exist if ’J'(X) = 0 vrhere i|)(X) is tslcen from equation (77) . 

Points of inflexion for the Figures 2.10-2,13 different 
values of B and Z are given below. 


FIG. 2.10 

B= .33,Xq=:0.1 

FIG. 2.1i 
B=0.5,X^=0.1 

FIG.' 2,12 

Scz 1 #0 • 1 

FIG. .2. 
B=2 ,0 

z 

Point of 

z 

Point of 

z 

Point of 

z Point 


Inflexion 


Inflexion 


Inflexion 

Infle 

.05 

,265 

.05 

.2794 

m 

o 

CJI 

,31203 

.05 .35E 

0.1 

.291 

0 .1 

.30 29 

0.1 

.3308 

0 .1 .369 

0.2 

.3494 

0.2 

.3568 

0.2 

.3755 

0.2 .404 


2 .12 G-EBERAL GLASS OF.MOIiBlB VaiH POMIS OF TOBEiXiaijr 


Besides the logistic and Smith* s models, we have given si 
more models of population growth based on essper imental laws 
connecting specific growth rate with the concentration of the 
rate- limiting substrate , Though these laws are different, 
yet they lead to similar population growth curves. This 
is not entirely unexpected behaviour. These laws are based 
on laboratory experiments with microorganism thriving on one 
rate-limiting substrate. Y/ith different microorganism and 
different substrate, the laws can be somewhat different. 

We now investigate mathematically general class of model 
with the same feature viz point of inflexion and limiting 
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We have got here a three-parameter family of curves 

Figures 2.10-2.13 show populations curves for different 
values of B and Z. Points of inflexion in these figures will 
exist if iJ»(X) = 0 where Tf»(X) is taken from eq.uation (7?). 

Points of inf lexical for the Figures 2.10—2.13 different 
values of B and Z are given below. 


FIG. 2.10 
B=,33,Xq=0.1 

FIG. 2.1i 
B=0*5 ,Xq»0* 1 

FIG.' 2.12 

B= 1 1 0 5 X ^ 

FIG. 2,13 
B=2,0,Xq=0.1 

7i 

Point of 

z 

Point of 

z 

Point of 

z P oint of 


Inflexion 


Inflexion 


Inflexion 

Inflexio 

*0 5 

.265 

.05 

,2794 

• 

o 

CJl 

.31203 

.05 .3559 

0.1 

.291 

o.i 

.3029 

0.1 

.3308 

0 .1 .3699 

0.2 

.3494 

0.2 

.3568 

0.2 

.3755 

0.2 .40438 


2 . 12 GBIJERAL GLASS OF. MOIELS Wl'fH PQMIS OF HTFIjEIIOT 


Besides the logistic and Smith’s models, we have given six 
more models of population growth based on e35)er imental laws 
connecting specific growth rate with the concentration of the 
rate— limiting substrate . Though these laws are different, 
yet they lead to similar population growth curves. This 
is not entirely unexpected behaviour. These laws are based 
pn laboratory experiments with microorganism thriving on cne 
rate-limiting substrate. With different microorganism and 
different substrate, the laws oan be somewhat different. 

We now investjgate mathematically general class of model 
with the same feature viz point of inflexion and limiting 


ANOTHER POPULATION MODEL 
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ANOTHER POPULATION MODEL 
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ANOTHER POPULATION MODEL 



Fig. 2.12 
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population sizes. Consider the model 

§ = g(H) (79) 

It will have a point of inflexion if 
2 

^ <g(H) + 8'(H) (H-JIell^=0 (80) 

is satisfied hy KH., where 0 < K < 1 
We consider some special cases 
(i) g(H) = AU+B 

Equation (80) gives 

Bh-A (KNg-^Ng) = 0 

or Mg(2IUl) + B « 0 (81) 

1 

B s= 0, K = ^ satisfies it, Ihis corresponds to the 
logistic law, 

Brom (81) 

I - OT- • ! = (1-8K)II^ (82) 

0 

-Q 

If ^ "then the point of inflexion occur before half 

ttie final population size is reached if B and A have the same 
sign and after this if B and A have opposite sign. 

If we know lig and K i.e, if we know the final population 
size and the position of the point of inflexion, then 

^ = A[:Mg(l-SK) 3 CMe 3 

gives the shiniest such model 


( 85 ) 



4V 

(ii) g(]!T) « + HSf + 0 (84) 

Equation. (80) gives 

iiE^e^+BKNe+O + (ZAKE^+B) (KlTe*^Te) == 0 (85) 

and we can choose any value of A,B,0 to satisfy (85) 

(tii) Similarly we can tafe g(n) is a pol^niomial of any 
degree or a rational function where denominator does not 
vanish between 0 and 1 

(iv) If g(H) = Ae®^^, we get (86) 

l+aEQ(I0-l) w 0 or a =i (87) 

Ihus 

St=Aexp[; ^ (88) 

(3U.K)N3 ^ ® 

gives a growth curve with limiting population size Eg and 
the point of Inflexion at KNe * 

¥e can thus get a large variety of models with points 
of inflexion and limiting population sizes. Even any data 
showing these features, can choose the model to give an 
optimum fit. 



CHAP'fEE - 3 


DII'I'EREl^GE EQTJATIOU MODELS IS BOOLOGI MD EPIilEMICS 

3.1 mmonjGTM 

Mathematical ecology and epldemology are among the most 
fa,scinating areas of study today. Gra/th and decay of 
populations, effects of pollution on population groY/th, prey- 
predator and host— paras ite models, conpetitian "between species, 
spread of epidemics and rumours a,ll provide fascinating 
applications of mathematics. Burghes [l9'75| demonstrated 
hoy/- population dynamics can "be used to motivate the study of 
differential equations at the undergraduate level. Kapur 
Q-978,197^ showed how problems of biomathematics can be used 
to motivate the study of partial differential equations, 
matrix theory, stability theory, stochastic processes, 
statistical mechanics, optimization techniques etc, because 
tiiese and other sophisticated tecimiques are needed to solve 
hundreds of challenging problems which remain unsolved in 
these fields. 

Differential equation models are accessible to caily those 
who have had a first course on differential and integral calculus 
and numerical analysis. Cm the other hand difference equation 
or recurrence relation models ere accessible immediately after 
a first course on elementary algebra and a few lectures on 


computer programming. 
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Both these types of models have "been used in the 
literature and "both can represent reality. In fact the 
validity of a model can only be tested by coinparing its 
predictions with actual observat icsns . Moreover differential 
equations are usually obtained by first writing difference 
equations and then using a limiting process. Bor numerical 
integration, vje have aga,in to transform the differential 
equations to difference equations. 

, In our discussions viie shall not even require the theory 
of difference equations. G-iven a recurrence relation 

and given a,b, we should be able to find successively 
u^jUpjUg,... • Similarly given the recurrence relations 

% = Vn-V ^4 

’ll = '>lVl+ Vn-1+ ^>4 (2) 

*11 = °lV.l+ °2VV °4 

and given a^, b^, c^ (i s= 1,2, 3,4), vie should be 

able to find out successively 

(u^jVijWf), 

with such siiaple equipment, we can get deep insight into 
some of the important problems of ecology, epidemology, 
economics, sociology, operations research etc. 
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3'.E DIE^PERENCE EQUATION MOIEliS 

In recent years there has been an increased interest 
in modelling ecosystems v/ith difference eq^uations, renewing 
an approach dating hack to ITicholson and Bailey [JL930 • 

This interest stems from the recognition that many species 
have important discrete con^onents in their life histories, 
such as synchronized breeding seasons and seasonal deaths. 

Difference equation models are more or less analogous 
to the more common differential equation models. Single 
species logistic models of various forms have been invest!.- 
gated by Maynard Smith []L96^, Pielou [jL9 6£j, Oook pL9?ll, 

May pL974,197^ and others. Parasite- host interactions, 

the original topic of the Hicholsan-Bailey investigation, 
have been 'Considered by Maynard Smith Q.968, 197^, Hassel 
and Varley [jL96^ and, with extensions of spatial variation 
by Allen p.97^ , Predator-prey systems have been studied by 
Maynard smith D96|], Maynard Smith and Slatkin C;i97g , 
Beddington et al Q.97^ , and others. Pretwell [1970 modelled 
species dynamics in seasonal environments by using difference 
equations Rolf Ifeller [1970 used difference equations for 
two prey-predator considering delayed population growth and 
starvation. Be examined how the different ass un^jt ions can 
alter the dynamic behaviour of prey-predator model, Adi Rave h 
and Uzi Ritte ]jL970 used difference equations for simulating 
frequency dependent changes in gene frequencies, Donald 1. DB 
Atigels Z19VQ coEpared birth and death model and a stochastic 
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difference ecLuation model, using data for a giant Canada- 
goose population. PHIL Diamond p.97Q used difference 
equation for estimating domains of non-global asymptotic 
stability and region of resilance for biological systems. 

Masao Halcumura pL97^ show that a certain limiting property 
of a system of linear homogeneous difference equations is 
characterized by a solution to a nonlinear system. Linear 
difference equations have been applied in estimation and 
prediction by Poole [l97^ and in age distribution models 
by Leslie 0-94^ . 

Ihe behaviour of difference equations in the neighbourhood 
of equilibria is known to be more complex than that of oorres-. 
ponding differential equation models used in describing 
continuous systems. Ihe choice between continuous and 
discrete population models from the point of view of stab ill tjj" 
has been discussed by a number of authors including Driessobe 
[jL 97£| , Duffin 10960, Lines 0-970, May 0975a, 19730 and 
Usher 0960 • It has been shown that the stability of 
discrete model represented by a system of difference equations, 
which is analogous to a continuous model represented by a 
system of differential equations, depends on the operator, used. 

May QsYsg shows that difference equations tend to be 
less stable than their differential equation twins, because 
the finite time lapse between generations of growth will have 
the destabilising effects associated with any 

A 62239 





5S 


interacting system. His discussion makes it explicit, 
clearly stability of the difference equaticns system 
ini? lies stability of the differential equation one, hut the 
converse is not necessarily true, Priesche pL974] using a 
particular operational rule has shown that the resulting 
analogous discrete system is as stable as the differential 
equation system. We shall obtain in the next chapter a 
class of operators for v/hich the difference equation system is mcc 
stable than the corresponding differential equation system, 

3 , 3 G-ROW'IH OP A PCPULA'IICasr 

We assume that the births and deaths in a unit interval 
of time are proportional to the population size at the 
beginning of this interval so that if X(t) is the population 
at time t, we get 

S(t4-1) - X(t) = BX(t) - DX(t) = aX(t) (4) 
or 

X(t+1) = (1+a) x(t) 

Given X(0), we easily find that 

X(l) = (l+a)X(0), X(2) = (1+a)^ X(0 . .,X(t)=( l+a)\(0) 

( 5 ) 

The graph of X(t)/X(o) is given in figure ‘3,1 for 
a — , 10 , ,15 and ,23 , 

It may be noted that the unit interval in the above 
discussion can, in practice, be chosen as small as we like 


GROWTH OF A POPULATION 

X(0) = 120 P X(0)=120<|, X(0)=t00 

a=0.25 I 0=0.15 1 a=0.1 
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by choosing "bhe time- sc ale appropriately, (Thus we may 
talce unit interval of time as a year, a month, a day, 
an hour, a minute, a second or even a microsecond depending 
on the population and the accuracy required, 

She smaller the unit interval chosen, the nearer will 
be the solution of the difference equation model to that of 
the corresponding differential equation model, but at the 
same time, the number of calculations required will also be 
larger • She unit time interval chosen thus may also depend 
on the computing device available which may be unaided 
mental calculations or logaritbmic tables or slide rule or 
pocket calculator as cn analogue cotiputer or a digital 
computer. 

In the present case, we may note that the corresponditig 
difference and differential equation models are 

X(t+1) - S(t) = aS(t), H =d]i(t) 
and their solutions are 

X(t) = (1+a)'^ X(0), X(t) = X(0) 

She two solutions are identical if 
e^ = ( 1+a) , 

3 GRaTSH 01' A PQPULASIdT WISH LBIISBD REBOURQSS 

Due to interval competition for limited resources, the 
increment in population size in a uni-b time interval is reduced 
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by an amount proportional to the square of the population 

size at the beginning of the time interval. If there are 

X(t) individuals, each conpetes for resources against 

everyone and as such the term due to competition depends 
2 

on X (t). Our model becomes 

X(t+1) - X(t) = aX(t) bX^(t) = X(t) [a-bX(t ) ]]] (6) 

Given X(o), a,b we can successively find X(l), X(2), 
X(3), . . .,X(t) , , . , , The graphs of X(t) for the following 
values of X(0), a,b^ are shovm in I'igure 3,2 


(i) 

X(0) = 30,0, 

II 

• 

o 

b « 

.0001 

(ii) 

X(0) = 10.0, 

3, — 

b = 

.0001 

(iii) 

1(0) - 5.0, 

Q, 5= #09 f 

b = 

.001 


It appears that the growth curves are all S-shaped and as 
t increases, X(t) approaches the value a/b. This can 
also be seen from (6) since when X(t) is very near to 
a /b, the increment in X(t) is very small. 

We can draw the graphs for other values of X(0), a 
and b. If X(0) is greater than a/b, the curve approaches 
the line X = a/b asyB|)totically from above . 

Sometimes the equation. (6) is replaced by 

X(t+1) = CX(t)/l+dX(t) (7) 

It can easily be verified that the graph of (7) is also 
S-shaped, As t 


\ &L 
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X('b) "*■ ( c<— l)/d. (8) 

In fact it can be verified by direct substitution that 
(7) has the exact solution 

X(t) = - .i:’ wWvi c •3*=- (9) 

dX^+ (c-l-dX^)c-*^ 

1 c»>W> c s i 

v/hich is a solution of a differential equation of the form 

^ = aX « bX^ (10) 

dt 

3 ,5 lUgLUMCB OF POLLUTION CM POPULA'JM GROWIH 

Pollution is a cumulative phenomena and depends on 
population growth. Pollution effect on population at 
time t depends in the populations at time t, at time t-1, 
at time t-2, and so on. Of course the pollution effect on 
pcpulation at time t of population at time t-I will 
decrease with I, V/e may tale it to be a decreasing function 
K(!I}) of I. 

We consider the model 

X(t+1)-X(t) = aX(t)-bX^(t)-cX(t)|X(o)K('t>)+X(l)K(t^l) 

+ ..+X(t)K(0);] (11) 

Given a,b,c, K(t), X(0),-we can find successively 1(1), 

4 - 

X(2), X(3),,.,, Pig, ,^.3, gives the graphs when K(t) = e" 
and 

Z(0) = 30.0, a = 0.04, b = 0,0001 
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and c takes the values 

c = 0.00006, 0.00001, 0.0002, 0.0003 

We note that as c increases, the population decreases • 

This is expected since increase of pollution should lead 
to decrease of population. We also find as t -► «>, 

Z(t) a finite limit which is verified to he 

( 12 ) 

h+c 1 .58c 

'We also note that the time taken to reach this value is 
almost independent of c. 

We may also take K(t) as a linear fimctim of t i.e. 

K(T) = 1-T/n, T < np K(T) « 0, T = n (13) 

so that 

x(t+.l)-x(t) = az(t)-hz^(t)~cx(t) |;;;x(t)+(l-. |) X(t~l)+ 

(1^ |)X(t-2)+..,. 

+ (1-|) X(t^)3 

(14) 

In this case, we should he given a,h,c, X(0), X(l),..., 
X(n.-1) and eq^uation (14) should he applied when 
ts=n, n+1, n+2,... Figure 3.4 illustrates the case when 
n = 4 and where a,h,c, X(o), X(l), X(2), X(3) have same 
values as in Figure ^3.3 • Bfere 

•* aa t - » (IB) 
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Eor the general value of n, 

t) ■* ■ . as t oo ( 16) 

Ih- ^ (n+l)c 

In model (11), we assumed that the polluticai before an 
arbitrary ' origin of time does not affect the future 
populaticm. This can only be 'true if there was no 
pollution before this time. V/e can also assume that it 
is only populations at n times viz. t, t-1, t-2,... 
t-.(n-l) which affect the population at time t so that 
we may consider the model 

X(t+1)--X(t) = aX(t)«bX^(t)-.cX(t)pL X(t)+e“‘'^X(t-l)+. . . 

+e“ (n-l)t(^jj:^)^ 

(17) 

In this case we should be given a,b,c, X(0), X( !),».., 
X(n«l) . It is easily seen that 


b+c(l-e“^)/(l-e"“^) 


(18) 


we ma.y note that in this case it is not possible to write 

4 

the analytical solution of (H), (l-l)* ( 1^) or the corres- 
ponding differential equation model 


= a3UbX -cX / K(T) X(t-T)dT 


(19) 


where I is the total time for which the effect, of 
0 

pollution is effective. I'or special values of K(l) 
(e.g. v/hen K(T)=:1) analytical soluticais however do exist. 
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Hhe difference equation recurrence-relation approach 
can however always he used to get insight into the influence 
of pollution. 

3 ,6 BIPIUMCE 01^ AGS SmUGOIURE 

For this Burghes [1.97^ obtained the time-lag differential 
equat ions 

11= e G^“ Yi;(t-T)3 x(t-a:)-Bx(t) (20) 

ror discussing its behaviour, he however used the difference 
equation 

X(t+1)-S(t) = ela- YX(t-I) 2 2(t-T)-BX(t) (21) 

and he drew the corresponding graph for 1 = 30, Xq= 5 and 
B = 0.05, ea = 0.25, 0 y= O.Ol, B = 0.05, 0a = 0 .025,6y=0.02 

B ss 0 .10, 0a = 0.25, 0Y= O.Ol? B = 0.15, 0a = 0.25, 0Y= 0.01 

( 22 ) 

To solve (21), we have either to specify the values of 
X(o), X( 1) , . . .,X(30) or we can assume as Burghes did that 
upto t = 30, only deaths tafees place so that (21) is 
replaced by 

X(t+l)"“X(t) — BX(t), ( t=0 j 1, 2 , » • . , 29) (2^) 

The exar^le illustrates hov'^ difference equations approach 
can be used to get insight into phenomena normally represented 
by delay differential equations. 
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3 .7 M AL'JHamSI^S LffiTHOD OF TAKING AGE STRUCTURE INTO 
AQO(gUNT 

We can divide the populaticai into three groups 

(i) Fafe-.jlepr oductive children (ii) Pre~repr oductive adults 
(iii) Post-reproductive old perscais 

Let X^(t), 

groups at tine t, then 

X^(t+1) — S2(t) = a^X2(t) — TOjX^(t) 

Xg(tfl) - XgCt) « agX^Ct) - "bgXgCt) (24) 

X3(iH-l) - X3(t) = agXgCt) - L^XgCt) 

In (24), the flcst term on the R,H,S. of the ficst eciuaticn 
represents number of births of children per unit time and 
this is proportional to the size of the adult group, and 
the second term represents deaths and migrations per unit 
time which are proportional to the size of the population 
of the children. In the same way agXg(t) and a 3 X 3 (t) 
represent migrations per unit time from the first group to 
the second and from the second to the third group and 
b X„(t) and b X, 3 .(t) represents deaths and migrations 
per unit time from the second and third groups respectively. 

Given X^(0), X2(0), 1^(0) we can proceed step by step 
to find X^(t), Xg(t), X 3 (t). 

We can write (24) in the form 


Xg(t), X3(t) be the populations of the three 



63 


, X(-b+l) = AX(t) 

where 


( 25 ) 



■'Xi('t) " 


-L-h^ 

Of 

0 

2(t) = 

Xg{t) 

A = 

^2 

l^bg 

0 


-V=>- 


_,0 

^3 



(26) 


SO that 

X(t) = a''^x(o) 


(27) 


The eigenvalues of A are 


^2 “ S" i 


Y(h^~h2)'+ 3 » 



l^hg 

(28) 


so that two of the eigenvalues are definitely less than 
unity and the third will he less than unity if ^^^2 ^ ^1^2* 
In this case the populations of the three groups will tend 


to zero as t -* ». If ^1^2* population of each 


group 

increases 

exp coient ially , 

. Irom (27) 



Xi(t) = 

t log 


t log Xg t 

log Xg 

fuOt: 


> 

°12® °13® 

■» A- 

X2(t) = 

t log 


t log Xg t 

log Xg 


°21® 


°22® °23® 


- n ; b i y 3 

3%(t) = 

t log 

Xf 

t log Xg t 

log X 



°32® + 



( 29 )*^ 

If we can observe X^(t), X 2 (t), X 2 (t) over a period of time, 
we can use the method of esponential peeling to estimate the 
values of the parameters a^jagya^jh^jh^jh^. 


64 


Figures 3. 5, 3, 6, 3. 7 give the growth curves for 2^(t) 



XgCt) 

when 

x^(o) = 

= loo, 

> Xg(0) = 

<* 

a 

o 

to 

(0) = 

100 


and in 










FI&. 

3.5_ 










(i) 


.05, 

#02 5 

br ■ 

,03, bg= 

•02 

^3= 

.02, 

bg= 

.05 





a^ag 







FIG. 

3.6 










(ii) 

a-=: 

1 

.04, 


^1= ' 

•03j 1^0=2 

/C> 

•04y 

^”3= 

.02, 

bg« 

.05 





a^ag 







PIG. 

3.7 










(iii) 

a^= 

.04, 

Q,p=: 1 0 2 j 

bl= ' 

•03j bg= 

•03, 


.02, 

bg = 

= .05 





^1^2 

^ ^1^2 







2 'S PREY-^REMfOR MD HOSS^iRASITB MODELS 

Let X(t), Y(t) he the populations of rats and cats 
at tirne t. Let us calculate the increments Z(t+1) - X(t) 
and Y(t+1) - Y(t) . In the absence of cats, population of 
rats increases and the increase is proportional to the size 
of this population, but due to the presence of cats, the rats 
population decreases and the decrease is assumed proportional 
to the product of the number of cats and rats. In the absence 
of rat , cats population decreases and the decrease is 
proportional to the number of cats and in the presence of rats, 
the cat population increases and the increase is proportional 
to the product of the number of rats and cats. Ihis conside- 
ration gives 



KO 



Fig. 3.6 
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(30) 


X(1H-1) - X(t) = a^X(t) - b^X(t) Y(t) 

Y(t+1) - Y(t) = -a2Y(t)+ bgXCt) Y(t) 

Given a^,ag,'b^,bp,X(o),X(o) we can calculate step by step 

X(t), Y(t) for any value of t, Figures 3 *8, 3 ,9, 3.10 
show the result of such a calculation for 

X(0) = 26, Y(0) =7,-fa^= ,01, ,02, b^-::.002, bg=,001 

XCol V ^ Vo,S 

'fhe equilibrium point for which X(t) = X(t+1), Y(t) = Y(t+-1) 


IS 


2^ = 20, Y3 = 5 

If the p op ulat ions start at these values, they stay at these 
value s . 


The populations shav periodic oscillat ians. In the 
phase space we do not get close trajectories as we get ia the 
case of the corresponding differential equation model. 


tl = ^1 ~ ■biX(t) Y(t) 

^ =-ag K-t) + tgX(t) Y(t) 


(31) 


However v^hen the unit interval is small, the trajectories are 
almost closed. Moreover observations of natural populations 
also do not give close trajectories. 

Here we have not considered self interaction. If this 
is taken into account then the equation (31) will become 

■X(t+1) - X(t) « a^X(t) - c^X^(t) ~ b^X(t) Y(t) 

2 (32) 

Y(t+1) - Y(t) -.agYCt)- CgY'^(t) + bgX(t) Y(t) 




Pred 



Time 


Fig. 3.9 
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Self interaotion is harmful te growth for both pop ulaticais • 
Behaviour in this situation for both prey and predater is 
shovm in i'igures 3.11, 3,12, 3,13 
where 

a^a .01 SLga .02 

,01 Oga .02 x(0) a 21.0 

b^a .002 bga ,001 x(0) a 6,0 

If there is a cover capable of protecting a limited number of 
prey, and of those prey unable to occupy this cover are killed, 
then all prey above some fixed number v/ould be killed. It 
may also iiappen that the predator is rare, and limited by some 
factor other than the abundance of this particular prey. 

Each predator might than take a fixed quantity of prey 
regardless of the abundances of the prey. In' this case a 
term X(0) is introduced in equation (31) which represents 
the number of prey those are not killed by predator. Model 
of this situation is given by 

X(t+1) - X(t) « aX(t) bl(t) [X(t) « 2(0)3 

(33) 

- XCt) = -oY('l;)+ dY(t) []C(t) -.2(0)3 

pqjulation growth of prey and predator is shown in Eigures 
3*1^, 3.15, 3.16 

where in sub figures of Eigures 3.14, 3.16, 3 *14, we take 


a 


.01 

2 a 


.01 

b 


.002 

b 


.004 

0 


.02 

0 


.03 

d 

ts 

.001 

d 

4- 

.002 
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3 8 , =: •04 

b = ,02 

0 = .03 

d t= ,002 


4 a = ,08 

b = ,009 

o *05 

d = .006 


X(0) = 21.0 and X(o) = 6.0 are coDnnon in. all sub figures. 
Both predator-prey populaticaa oscillate with time, The 
behaviour of prey with respect to predator is shown in 
Bigure 3.16 . It shows that as time increases they tend 
to equolibrium value , 


3.9 IHTBRAGTIOH BBTWEM PEEBATCRS MB PREY 

Le-fc are the number- of prey and are adult 

and young predators. In each interval, each predator must 

find and eat one prey. If it fails to find a prey it dies. 

a, 

If it succeeds it does not seek fer second prey until the 
next interval. The model of this type is 

^t+1 ” ^t ^ 

^t+1 ^ ^t+1" ^t+1 

Here a, and p are the effective search areas of adult and 
young respectively. 


If we plot a graph between prey population and time, we 
see that prey population decreases exponentially and after 
some time it will become constant. The behaviour of adult 



50 


Time — 
Fig. 3.11 


100 


150 
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and young predators is similar. Both populations decrease 
exponentially and will die out ultimately. Behaviour of 
and is shown in 

Fig * 3.17 F ig . 3 .18 Fig • 3 . 19 

She values of constants and initial values are 

a = ,006 Z(o) « 400 

Z(0) = 40 

P = '.004 Y(0) = 60 

3 .10 IiESlIB»s EOTJAIIOES 

Xtt -rO aXft:;- hX^’ ~ 

» - 2 / 

ya-\\) - fY^/XCtJ 

This is the prey-predator eq^uation with some modification. 
First eq.uation is of prey and is the same as Yoltera*s 
eq.uati(m with damping. Modification occurs in the second 
equation. If prey per predator i.e, x/Y is large then 

V 

predator increases exponentially if X/Y = f/e then the 
predator is a.t equilibrium. If X/Y falls below f/e 
the predator will die out. 

Behaviour of this equation for 


a — .0 1 e 

= .002 

O 

• 

1! 



Y(0) = 3.0 

b =: .02 e 

= .02 ' 

O 

11 

• 

o 

shown in Fig . 

3 ,20 ’ and 

Fig • “ 3 .21 , 


jjf- Sat S' ftW'ywJuy 
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Prey population decreases e 2 q) onent ially and finally becomes 
extinct. Predator first increases upto certain value, 
after that it will decrease sharply and then again start 
increasing slowly, 

3.11 CCMPEPITIOH mWEM WO SPECIES 

In this, each species increases in the absence of the 
other by an amount proportional to its own population, but in’ 
the presence of the other, it decreases by an amount propor- 
tional to the product of the two so that our model is 

X(ti-l) - X(t) = a^x(t) - b.X(t) Y(t) 

(36) 

Y(tfl) - Y(t) = a^YCt) - bgX(t) Y(t) 

* 

Figures 3.22, 3.23, 3,24' give the graphs for the case when 
a^ = '.Ol, ag = .02, b^ = .002, bg = .001 so that the eq.uilibrium 
point is X_ = 2o, Y = 5. 

The survival of a species depends critically on initial 
size of the two populations. When X(o) = 25, Y(0) — V or 
X(0) =5 15, Y(0) = 3.5, the first species tends to disappear 
and when X(0) = 25, Y(o) = 6 or X(0) = 15, Y(0) = 3, the 
second species tend .to disappear , When X(0) = 25, a change 
of Y(o) from 6 to 7 makes the seccnd species survive and the 
first species to disappear. A similar effect is seen when 
X(0) « 15 and Y(0) changes from 3 to 3.5. 
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If oat of the two species, only one is of the desirahle 
type, we can ensure its survival and the vanishing of the 
other hy just adjusting the initial populations. 

If self interaction is also taken into account then 
the above model represented by the equation (36) v/ill become 

X(t 4 -l)-][(t) « a.jX(t) - b^X^(t) - c^X(t) I(t) 

2 (37) 

Y('b+1)-T(t) « agY(t) - bgY'^(t) - C 2 X(t) Y(t) 

both species will decrease by an amount proportional to 
the product of its own population. 

Population growth of first coii^eting species is shown 
in figure S,25 and of second in Pigure '3.26 while the 
combined behaviour of both species with respect to time is 
shown in Pigure 3.27' • 


In sub figure 

i 


a^ = .01 

ag = .02 

x( 0 ) = 15,0 

b^ = .02 

bg = .03 

Y(0) 3.0 

=5 .002 

C^ a .001 



In this case second species will grow while first species 

will die out • 

In sub figure 2, 


= .01 

a^ = .02 

x(0) = 21,0 

b'^ = .005 

bg = .008 

Y(0) = 6.0 

c^ = .002 

0 g .0 0 1 


In this case both species will die 

out. 
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3 ,12 GOUPETITIOH MOUG IHPJaS SPECIES 

In, this model, each species grcws in the absence of 
others, hut mutual contacts are harmful since they are 
competing for the same resources. Our model can also 
represent tliree tribes at war with each other.. The 
model is 

X(-tH-l)-2(t) a^X(t)-agI(t)Y(t) - agX(t)Z(t) 

Y(-b+l)-.Y(t) = b^Y(t)-b2Y(t)2:(t) bgY(t)Z(t) (38) 

Z(tfl)-Z(t) = c^Z(t)-C2Z(t)X(t) - CgZ(t)Y(t) 

The graphs below corresponding to t''ie following values 

— *01 f ^2” ,0015 a^ ~ ,002, b^ — ,002, bg = ,002 

"bg = .003, c^=o,03 , c^ =0.03, Cg = 0.004 

are shown in figures ^3.28-3.32 
The equilibrium point is 

Xg = 0.6, Yg = 3,4, Zg = 3,6 

sub figure. • corresponds to X(0) = 5,0, Y(0)=7,0, Z(0)=7,0 

sub figure. correspondsto X(0) = 7,0, Y(0)=:l0.0 , Z(0)=10»0 

It appears that for case (i), the first two species 
disappear and the third species alone survives, while in 
case (ii) with different initial values, the first and third 
species tend to disappear and the second species alone 
survives. This again . illustrates how initial population 
sizes oan play a critical 3 ||Dle in the survival of the species. 
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In. tiiis model each species gxov/s in the absence of other, 
hut mutual contacts with rem;un±ng two are useful for first 
species* For second species, mutual contact with third is 
useful hut harmful to first. Mutual contacts are harmful 
to the third species with remaining two. The model is 

X(t+l)— X(t) = AX( t )-hBX( t )Y( t ) + cX(t)Z(t) 

y(t+l)-Y(t) = DY(t)+EY(t)Z(t) - FX(t)Y(t) (39) 

Z(t+1)-Z(t) = GZ(t)-HY(t)Z(t) -EX(t)Z(t) 

The graphs corr e sp ending to the folio-wing values 

A « .01 B « .001 0 « .002 D a .04 E = .002 

E « .003 G = .03 Ha .003 K = .008 

are shown in Eigures 3.33-3.30 

Suh Figure 1 corresponds to X(0)=35.0, Y(o) = 2o.O 2(0 ) = 35.0 

Suh Eigure 2 correspcnds to X(0)=25.0, Y(0) = 9.0, Z(o) = 25.0 

It appears for both cases, the first species only survives 
while second and third v/ill disappear. 

EPTUEMIG MOEBL WITHOUT PiEM OVAli 

Let X(t) denote the number of susceptible and Y(t) 
infectives at time t. In the time interval (t,t+l), X(t) 
decreases by px(t) Y(t) and Y(t) increases by the same amount 

so that 

(^) 


X(t+1)-X(t) = -pX(t)Y(t) 
Y(t+1)-Y(t) == ^X(t)Y(t) 






0 7 






«r-- 


Third species Z(t) 



:OMPLT!TiON AMONG THREE 
:ONTACT lO HARMFUL for" 
-OR f'-'THEIR 


FLCiE'-' 
1 GPE 


WHERE MUTUAL 
JL-- AND ■..UCFUL 



First species X(t) — > 
Fig, 3.38 
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Initially let there he n sus cep titles and 1 infective, 
so that 

X(0) = n, Y(o) = 1 

Prom (4o) 


X(t+l)+Y(-tfl) = Z(t)+Y(t) = l(t-l)+Y(t-.l)=... X(0)+Y(0)=n+1 

( 41 ) 

so that the trajectories in the X— Y plane are straight lines* 
Prom (40) and (41) 

X(t+l)=X(t)-px(t)[n+l^X(t2= X(t) [l-pn-0+ px^(t), 

(42) * 

so that 

^ X(l) = n-np, X(2) = n(l-.p) (l^p-np^) (43) 

and we can successively find X(t) for any given value of t. 

The graphs of X(t) and Y(t) are given in Pigures 3,39 and 

3,40 for n = 10, p = .05, ,03, ,02, ,01. 

It appears that the number of sucep titles falls rapidly 
and the number of infect ives rises rapidly and the greater 
the value of p, the greater is the rate of fall or rise, 

3 .14 LPIMIG MOIEIi WITH REMOYAl 

In this case (4o) becomes 


X(t+1)-X(t) = -pX(t)Y(t) 
Y(t+1)-Y(t) = pX(t)Y(t),, YY(t) 

so that 

It jvirt-ii) “7 Ot 


(44)Y-;(t 
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X(t^-l)+Y(i>^-l) = X(t) + Y(t) - YY(t) (45) 

Yhe graphs of X(t) and Y(t) are given below for 

n = lO, p = .01, .05j .1, Y= 0.1, 0.2 

figure 3.41 Figure 5,4^ 

% 

It appears that the number of susceptible fall and the 
rate of fall increases with both p and y . \Yten p is 
large, it appears that the number of infective rises 
quickly and then falls. 

3.15 G- aOV/IH OF PQPTJlAflOH WITH IME-IDELAY 

Growth of a population with time-delay effect when 
abundant of food is available is studied. Ihe model is 

X(t+1)-X(t) = aX(t-30) (46) 

Behaviour of the population growth with respect to tinas is 
shewn in fig, [3.43 for a == ,0l and in Fig, 3,44 for 
a — .04 . 

For sinjtlicity v/e consider that populaticai remain 
canstant i,e, X(t) = lo upto t = 30. After t = 30 it will 
grow e 2 !panentially, Ihe growth of the population is similar 
as in the model where time— delay is not considered, 

3.16 GY07/TH OF PaPUMIION WIIH IBIB-BEMY Y/IIHOUI 

gtcyroiiiy Pbos 

Growth of a population when sufficient food is not available 
and tte effect of past is considered, is studied, Ihe model is 
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EPIDEMIC MODEL WITH REMOVAL 



Fig 3.41 



CjKUWTH of population with time- delay 104 
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X(t+1)-X(t) = -.aX(t) ~ M(t-30) (47) 

Due to shortage of food populaticai will decrease# 

Behaviour of the population with respect to time is 
shown in Pig. 3,45 , 

In suh Pig. 1 effect of shortage of food per unit 
population is the same as the time-delay effect 
on unit population i.e, 

a = .001 and h =* .001 

In rest of ' suh figures effect of the shortage of food 
per unit population is more than the effect of past on unit 
population 

i.e. in Suh Pig. 2 a *= .01 and h ^ .008 
Suh Pig. 3 a = .01 and h « >005 

Por simplicity we ccoisider that population remains constant 
i.e. X(t) = 10 upto t = 30. After t = 30 it will 

decrease and vanish* 

3 , 17 pESy -PRP.'nATOR MQDEli WITH TME-DEMy ^ 

It is possible that the reproductive rate depends, not 
only on the population density at the present time, 
the population density in the past. Ihe reproductive. of 
a predator species will depend on the prey, which may in turn 

depend on hovf muoh of the prey was eaten hy predator an the 
last years. To gain some idea of the effect of such a 
delay in the effects of population density a model is given hy 

L 



106 



Fig. 3.45 
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X(-tH-l)-X(t) = aX(t) ■bX(t) I(t) 
l(t+l)_Y(t) = -.cl(t)+ dX(-U20) Y(t-20) 


(48) 


where X is the prey populatissn and Y is the predator 
p opulaticai • 


Behaviour of the prey, predator population with respect 
to time is shown inl'ig. ,2.4:6 and Pig. 3,47 respectively. 
Pig. 3.46 shows mutual hehaviour of prey-predator 

in Suh Pig. 1 a=.Ol, 13=. 002, Cs.02, d=»001 
Suh Pig, 2 a=.0l, h=.004, c=:.03, d^.002 

In hoth cases prey populaticsn ?/ill in-ore ase vfhile .the p-^©4rertor 
population die cSut. ^ the time increases. ^ 




T ime — » 


Fig 3.46 



PRLY - PrYLDATOR MODEL WITH TIME DELAY 



Time > 

Fig. 3.47 




OHAPIER « 4 


RBIiil2Il/E SIABILin OP PHiTEtffiNGE MD MPPERElTTIAl 
EQUATION SYSTUIS 

4.1 I NTRODUCTICN 

The choice "between continuous and discrete population 
models from the point of view of stability has been discussed 
by a nuxflber of authors including Dries sc he 1^19 V‘|] » Duff in 
[ 1069 ], Innes 0-974], May 0-9 7|1 and Usher . In 

particular May has considered the continuous model 

represented by 

^ = P.[|,(t), N„(t)...N (t)0 j 3 = 1 , 2 , ...,m (1) 

dt 3 ^ ^ ® 

and the fmalogous discrete model 

N .(t+l)-.N.(t) = E.[N.,(t), N2(t)...N^(t)3l 3=I»2,...,m 
3 j J 

( 2 ) 

where N^(t), N 2 (t) , . . .,Njt) are the populations of the 
m species and P^*s are non-linear interaction functions. 

Por defining the analogous discrete system, May has replaced 
the derivative by a forward difference operator with step 
size h > 0 so that 

N^(t) ^ ^ (T-1) 

where 

TN.(t) = N^ (tfh) 

D 3 

and time has been normalised so that h — 1 


( 4 ) 
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In the two models, the biologloal featurao such as 
teophio etruotures, birth and death rates, competition and 
prey-predator interaotions are identloal. The equilibrium 
pcpulaticns are also the same. The main difference is 
that while for (E), growth takes place in discrete steps, 
an (1) it is continuous process. May [igyg has stated that! 

* 'It is vifidely understood that difference equations 
tend to be less stable than their differential 
equation twins, because the finite time lapse 
betv/een generation of grov/th will have the destabi- 
lizing effects associated with any time lag in an 
interacting system. Our discussion maJces it explicit, 
clearly stability of the difference equations system 
implies stability of the differential equation one, 
but the converse is not necessarily true*’ , 

However Driessche QsVfj , using the operational rule 

(®-l) ?: I (3h-l) 5^-- (5) 

has shown that the resulting analogous discrete system is 
as stable as the differential equation system. 

Here, we examine the operational rule 

h 

= pfi 3/ 

v\fhere (>( -1) is a parameter, to compare the stability of 

the continuous and discrete models. It' is obvious that when 
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p = 1, (6) reduces to (5) and when p = 0, (6) reduces to (3), 
so that the results of Drlessche (1974) and May (1973) should 
follow as particular cases of our results and v/e may espect 
to get additional insight hy using other values of p . 


4.2 MaLYTIOAL hISGUSSIOI 


Let n* >0, (3 = l,2,,..,m) he the equilibrium 
p op ulat ion value s , whe re 

Sj (N* , kJ N*) = 0 (7) 

3)0 discuss the stability of the equilibrium pop ulat ton, we 
perturb these by' small disturbances 

IT.(t) = N? QL + V.(t) ] ( 3 = 1 , 2 , .. .,m) (8) 

D 3 J 

Substituting in (1), neglecting squares, products and higher 

pov/ers of V.(t) and using (7), we get 
^ 3 


- S (d= 1 ? 2 , . . .,m) 

fc=l ^ 


where 


a 

- j,* ^ 911^^ 

D 


(9) 


( 10 ) 


and the partial derivative are evaluated at the equilibrium 
population values. Substituting 

, V -n (1^) 

V.(t) = B. e 

In (9), we find that this can be a solution provided 

, ( 12 ) 

liU XI) = 0 
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where a is the m x a matrix The system (1) will 

be staole if for all eigenvalues of A, we have 

B-e ^ < 0 (13) 

Substitutitig (8) in (6), we get 

, V h dV. 

V.(t^h) - V.(t) (1+ PT) ^ (14) 

using (9) 

Tj('tH-h)_V.(t) = 2 (16) 

Substituting 

V^(t) = Oj (16) 

v/e find that this will be a solution provided it satisfies 
the etiuation 

The analogous discrete model vrill be stable if 

I'J' I < 1 (18) 


4.5 G-EavIEiTrllOAl IIO!EHPKSTATI(|[S 


CotEparing (12) and (17) ^'e get 


, ?+l t^-1 


(19) 


or 


, h 1+P+ ^ h 


(20) 


‘ibB olroular regim 1 ♦ 1 < 1 correspmda to tto reglco 
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1+P+ X h| < 1+p— p X]i| 


( l+p+Xh) ( 1+P+ Xh) < (l+p~pXh)(l+p— ph x) 


XX h^(l-p^) + Ii(i+p)^ (x+ X) < 0 


(xW) (i-p^) + ^ (l+p)^ < 0 
If 1 pj < 1, (21) givBs 




( 21 ) 


( 22 ) 


which is t]ae interior of a circle with centre (- ^ ,0), 

radius r an<i which passes through the origin, 

h 1- P 


radius 


(ii) If jpj > 1, (21) gives 


fY ^ P+l ') 2 S ( P t 1 a ^ 


(23) 


1 P + 1 


v^hich is the exterior of the circle vrith centre (j- » O) 

radius — £> 3 ^d which passes through the origin, 

h p-1 

(iii) If p=^ 1, (21) gives 


X < 0 


(24) 


which is the left-hand half of the X plane. 


We show in I'ig. 4.1, for different values of p , the 
regions corresponding to j'/'j < 1* 

We get the following possibilities, 

(i) When p = 0, we get the interior of the circle of 
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1 

centre (- ^ ,o) and radius^ » This corresponds to the 
Case discussed by May 0.97;^ , 

(ii) As P increases, v/e get interiors of larger cond larger 
circles with centres on negative x-axis and all passing tlnrough 
the origin, 

(iii) ',/hen Pel, wo get the v^hole of the left hand plane. 

This corresnands to the case discussed by Driessche [197^ . 

(iv) v/hen -1 < p < o, we get interiors of circles of raxlius 
less than - v/ith centres on negative x-axis and passing 
through the origin. As p -1 + o, we get just the point 
circl:: at the origin, 

(v) * When p is slightly greater than unity, we get the 
exterior of a large circle with centre on tl^ positive x-'’xls 
and passing through the origin. This exterior includes tlie 
entire left hand half of the x-y plane and part of the right 
hand half of the x-y plane . 

(vi) As P increases beyond unity, we get exteriors of 
smaller and smaller circles with centres on positive s^axis 
and passing through the origin, 

(vii) A P -*•«>, we get the exterior of a circle of centre 

1 1 
(■^ , 0) and radius ^ , 

(viii) V/hen p in negative and large, we get the exterior of a 

-t ... 

circle of radius Isafg&r' than , centre on the positive xj-axis 
and passing through the origin, 

(ix) As p - -1-0, we get the circle of zero radius at the 
origin, ’ 
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(30 As h 0, both the regions given by (22) and (23) 
^procich the region X < o and in the limit both system 
tend to be Gq.uistable • 'JJhis is of course expected. 

4 .4 OOlIPARISOtT OP THE STABIIIIY OP CON'flNUOTJS MS 
pTsaRE^-iroSPijliS^ ■ — — — — 

Oa^se_^I -1 < P < 1 

In this Case if the difference equation system is 
stable, I if; I <1 and the corresponding point in the X -plane 
lies v^ithin the ciccle v/ith centre (- ^ , o) and radius 

K ’^"’p ^ shaded region I so that X < 0 and the 

differential equation system is also stable. 

On tbe other hand if the differential equation system 
is stable, IL X < 0 and the corresponding points lie in tlie 
lef't-hand half of the X -plane, but this point can be outside 
shaded region I and |i|<l can be greater than unity so that the 
difference equation can be unstable. 

Thus in this case the stability of the difference equation 
system implies the stability of the differential equation 
system but the converse is not necessarily true. 

Case II P - 1 

In this case jif;! <1 => R 3 X < o' and R^X < 0 ==>l.f;l<l 
so that the two systems are both stable or both unstable. Ihis 
is the case discussed by hriessche p.97|] . 
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0ase_ III 1 p| > 1 

In this case the region j'Pj < l corresponds to the 
region III in the x plane. When Ul < 1, R x may he 
positive or negative. Thus R^ X < o => |i|»| < 1, hut 

I'^l ^ ^ ^ ^ 0 so that the stability of the differential 

eq_u.aijion system, implies the stability of the difference eq^ua-tion 
system, but the converse is not necessarily true, 

4 .5 DISGUSSION' 

we have shovn that the relative stability of the 
difference and differential equations models depends on the 
operator used and we have demonstrated an operator for which 
difference equation model is actually more stable. 

When P < 1, difference equation model is less stable, 
but by making h smaller and smaller, we can make region of 
stability larger and larger till it coincides with the region 
of stability for the differential equation. 

'./ben P > 1, differential equaticn is less stable, but again 
.by making h smaller and smaller, we can make its region of 
stability larger and larger till it coincides v/ith the region 
of stability of the difference equation, 

>’/hen p = 1, the region of stability is the same for all 
small values of h , 

This discussion is illustrated with an age-structured 
population model with three age-groups* 
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4.6 M A&E~STRUO[i]X]IlED PCP ULATIOIT MQIEl 

\/e divide the population into three groups, the 
pre— reproductive , the productive and the post— reproductive . 
let d^jdgjdg he the dea.th rates in the three groups and 
bg be the birth rate in the reproductive group. Also let 
m^jUp be the rates at which individual from first and second 
groups migrate to second and third groups respectively, then 
the basic differential equaticn for the model are 


or 


where 


dXg 

dt“ = "^1^1 " 

dX„ 

— « mgXg - dgXg 


(25) 


ax 

dt 


= AX 


(26) 



X^(t) 


*“( d^+^-j^) 

bg 0 

x(t) = 

Xg(t) 

x^i^) 

5 A s= 

mi 

0 

-(dg+nig) 0 

mg -dg 


The eigenvalues of A are given by 

4( d^+m^) ( V“2)‘^^^2°^l ^ ^ 

= i {-(d^+m.+dg+^Pg) + C (^i^“l"^2‘'“2^ 4bgm^3 ^ > 


(27) 

f2R’) 
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All roots are real. iVo roots are negative and the third 
is ivs or positivo according as 

(11+ m^) (dgH- mg) I bgm^ (29) 

Thus the equilibrium position (0,0,0) for the model (25) 
would be stable if 

(di+ m^) (dg+ mg) > bgm^ (30) 

and it vo\xld, be unstable if 

(d^+ mp (dg+ mg) < bgm^ ( 31 ) 

For this model, all the eigenvalues are real and therefore w© 
discuss in the next section, the relative stability of 
differential and difference equation model for the special 
oa,se when the eigenvalues of the differential equation model 
are real, . 


4 .7 xiEh/iLTIVE STABILITY 01? in'TStMTIxiL ML LDTEEiENCE 

mrmmrtsmmimr — - — 


prom (20), K= Xh 


Oa^e^ 

J. 

^ > 0, Xh — K > 0 


Ocaisider the curve 



li-x+K P, 

y “ 1+x-Kx * 

When 

jL 

= 0 , y = i+K > 1 

when 

X 

Y > 1+K 

when 

X 



(32) 


(33) 
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The curve has tv^o asynptates 


X = -, -i- , X ~ -L 


when 

11 

0 

11 

1 


when 

Y = 1, X = 


when 

V __ *1 '\r 2-l"K 


If 

- ^ ^ < -1 . 

then -1 < Y < 1 

If 

“ izra < *’ < -1. 

then -1 < < 1 or 

If 

P >-1 or P < - 

' 1^1 > 1 


(34) 
L 

(35) 

(36) 


The curve sketched is shown in kig* 4.3 
^ < 0 , ^ h = _k, k > 0 


\ie consider the curve 


l+x~k 

1+x+kx 


0 < k < 1 


(37) 


proceeding as above, v?e get the curve sketched in Figure 4,4 
When *«1 < X < — , |Y| >1 


When -1 < p < - , I'l'l > 1 

when p<-l or p>-. > I ’f' 1 < 1 • 

Thus we get the situation represented in the Figure 4,5 
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1 > 0 

2+Xn _.l 0 


1 < 0 LiL> 1 

-1 ^+Xli 0 

-SITET 

Figure 4 ,5 i 

'ihus \'/hen X is real, we get the following results 

X 

(i) 'ihen P < min (-1, ^ or p > max (-1>- > 

either both models are stable or both are unstable, 

(ii) 'When p lies between and ~ , one modfel is 

sta.ble and the other is unstable . 


'Dhus when X is real, we find thau when P lies betv^een 
,-1 and - , either X > o and < 1 or X< o 

and |i|)l >1, 


4 ,8 EXBjgliB OF AGS-SFRU GlT mED PQPTTIA'IIOH MOIEL 


Case 1 bgffi^ > (d^+m^) 

There are three eigen values — X^, “^ 2>^3 where 


X 


ijXgjXg 

If 


> 0 . 

p is taJcen between 



2+Xgh 

o 


and < 1, 


In this case the differential eauation model is 
unstable and the difference equation model is stable. For 
other values of p , both models are unstable , 
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tfere the three eigenvalues are ^x «x 

X ti* o 

If P is talcen in the intersectioa of the interval 

2 ->^-h 2-X h 2~x h 

s+qn » (- (" z+igL^ 

then Ill'll > 1 , l-j-gl > 1, I'Pjl > 1, 

In this case the diffrential equation is stable and the 
difference equation model is unstable. Por other values 
of p, both models would be stable, 

Ihe smaller h is, the smaller is tbe interval for 
exceptional values of p for which both models behave 
differently. When h -* 0, both models are stable or unstable. 



GilAPTER - 5 


BIFURGATIOH THEORY i'QR TwO PCXPUIATIOH MOHSLS 


G.l ] MiOj) UOTIOH 

Bifurcation theory is a study of the branching of 
solutions of nonlinear differential eq_uation. It is of 
particular interest in bifur cat icaa theory to study how the 
solucions and bheir multiplicity changes as parameters 
vary* Thus at a bifurcation point there is a transition 
in solution multiplicity. Application of Bifurcation 
theory are in elastic and hydrodynamic stability. In 
elo.stic stability theory, bifurcation is called budding. 

In hydrodynamic stability, bifurcation points are called 
transition or critical points. 

Bifurcation theory, loiown as branching theory as a 
mathematical discipline, had its origin towards the end of 
19th century, in the work of Poincare (1892-99) who made 
use of nonlinear perturbation theory in celestial mechanics 
and who was followed by Liapunov 0-900 and Schmidt p.91(2 
v\rho, in turn worked in the area of nonlinear integral equations 
and then, in la,te ti\renties by Bride brandt and Graves |_1920 
provided an essential step in the development of bifurcation 
the ory , 

Important developments v^rere made by the Soviet school of 
Krasnoselsku [JL964| and Vainberg lls60 and significant gravth 
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was experienced in the sixties. Ihe basic mathematical 
ideas are contained in the booh of Pimbley p.96g while 
numerous physical applicaticms can be found in the booh 
edited by Keller and Antaman [196^ • Bifurcation theory 
made significant progress aad has find application in 
various estaolished fields and to many new and emerging 
sciences. Bifurcation theory has been widely acclaimed and 
has found fruitful applications in theoretical biology. 

J .P.G. Auchmuty and G* Kicolis pL975] using bifurcation 
theory constructed new steady state solutions of a nonlinear 
ne'bvYorlc involving chemical reaction and diffusion. Same 
anthers in [jL97^ described oscillator solutions of a system 
of re aotioru- diffusion equations. Bere solutions exist for 
a wide range of values of the parameters and may be calculated 
by using the technique of bifurcation theory, M. Herschhowitz- 
Kaufman GLgvC made a detailed comparison for different 
boundary conditions between the analytical solution of the 
Kinetic equations obtained by bifurcation theory and the 
results of computer simulations. \7.H, Pleming's [197C1 models 
on geneticsj deal with selection migration and his result 
ccxnoern the exis'tsnce and stability properties and bifurcation 
phenomena for equulibrium solutions* J. Guckenheimer and 
G. Oste-i^ s’^o.died dynamics of density-dependent population 

models. They illustrated the general theory on a density- 
dependent Leslie model with two age classes. Ihe pattern of 
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bifuroatim way from tto aqullltelum point was invaatigated 
rieh IfeWian [jays: uaed too ^thode to study stability 
criterion for the Hopf bifuroation of dlffeomorphistas cm 
He lormula was then applied to the study of a density 
dependent Leslie model of tv/o age classes. George H# 


PimLley p-G?!] in a pr 6 dator-pre 3 ’' model allows that for certain 
parametric values there exist a closed separatrix in the phase 
plane and gave conditions under which this closed sepratrix 
is suable or unstable. Then he used bifurcation theorem of 
Andronov and Leontovich to shov/ how larger amplitude periodic 
solutions branch from this dosed sepratrix. J.F, Perez pL97^ 
disciiosed population oscillations in isolated populations of 
flies . Ife analyzed the solutions in terms of the parameter 
range . Ife found a condition which shows that a pair of root 
crosses tiie imaginary axis a single time. Each time a pair 
of roots crosses the imaginary axis; qualitative changes are 
expected in the solutions where this occurs. -Ibis is called 
Hcpf bifurcation. J.II. Gushing p.9VCl considered a general 
integr odifferential system which described a predator -prey 
interaction subject to delay effects. Ife shows that there 
are critical values of birth and death of prey and predator 


at which stable periodic solution bifurcate from tlK equili*. 
brium. He does not consider the stability of the bifurcating 
periodic solutions . Sle problem of determining theifebility of anon.- 
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constant periodic solution is in general quite difficult for 
an-cegrodifferential systems (and even for ordiuany differential 
systems). In this regard MacDonald [197^ considered the 
hifurcabion of periodic solutions for a n=2 predator— prey 
model with delays. Mac Donald’s approach is to convert the 
integr 0 differential system to a larger differential system 
and then apply Hopf hifuroaticn techniques. These results 
of Mac Donald for predator— prey model are analytically 
obtained stability results for the "bifurcating periodic 
solutions of integrodifferential systems. J.M. Gushing 
[1.977a] asserts the existence of nonconstant periodic solution 
which bifurcate from equilibrium as certain parameters pass 
tlicough critical values. J.M. Gushing [l977'Q modifed predator- 
prey dynamics represented by differential equations by the 
assumption that the coefficients are periodic function of 
times. By using bifurcation theory he has shown that this 
system has a periodic solution. Numerical results which 
clearly show this bifurcation phenomena are briefly discussed. 

In tlie present chapter, a set of conditions for a pair of 
eigenvalues of a stability problem matrix to be purely 
imaginary are obtained for a general physical system* These 
conditions are then applied to discuss the existence of 
bifurcation points for (i) time-delay prey-predator model and 
(ii) gro\>?th of a single population under the effect of poUu'cron. 
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5.S or POBITS OP BIPURCAIDIOIT A GENERAL 

pMlT m TIEW — ' — — 

Consider a physical system goYerned hy the systern of 
dix f e rent irl e quat ions 

dX. 

dt" (i=l,2,...,n) (1) 

where ^2’ * “ ’ certain parameters • 

Its equilibrium points are given hy 

^i^^l*^2^ * * * ^'^^^ ^1*^2^ ** *’^m^ ~ ^ (i=lj2j , . » jU) (a^) 

Lot (i^jXg, . he an equilibrium point . 


Let 


X.. 


,• + h. 


(irslj 2y * • • jU) 


(3) 


'i r 1 

Substituting in (1), using (2) and neglecting squares, products 
and higher povrer we get 


5 

dt 



n 


D=l'^3 


u. = 2 a..u. (say) (^=1,2, . , .,n) (4) 

5 3=1 3 


Substituting 

== A.e^^ (i=l,2,...,3i) 

i 1 

v/e find that (5) will give a solution of (4) rf 

liUXll = 0 

v^hore A is the matrix 


(5) 

( 6 ) 
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^11 

^12 

In 

^21 

^■22 

^2n 

• • • 

« • • 


« • • 


• •••« 

a ^ 
nl 

®'n2 


^ • m 


Let equ, (6) Le written as 

K2x“-®+..wK„= 0 (8) 

v/hore K^jkg, . ,,,K^ are functions of tlie parameters C^, Cg, . . . , 0^. 

Lor discussing bifurcation, we are interested in finding 
whether (7) or (8) can have purely imaginary roots of the 
form + iai . Ihe interest arises because the nature of the 

MM 

equilibrium point changes when parameter values change 
slightly in this case . For the two dimens icaaal case if the 
roots are y + i«») , we get an unstable focus if 7 > 0, a 
centre if 7 = 0 and a stable node if 7 < 0 . 7 is a 

function of the parameters and for some values of the parameters, 
7 can be zero, A slight change in the values of the parameters 
would change the nature of the equilibrium. 

If + i <0 are roots of (8), get 

...) t = 0 (9) 

Equating real and irflag inary parts seperately to zero, we get 
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n ,, n-4 , (»-l)/2 

u) +^4“ -...+(-1) ^11^1“= 0 


IC^»"-l-K3«i“-®+K5»“-6_..^(_l)(n-l)/2K^ ^ 0 
if 21 is odd aiid 


n 


( 10 ) 


w^-Kg ( - 0 

I^__^a)=: 0 


( 11 ) 


is n is even. 


For (8) io have a pair of purely imaginary roots, equations (lo) 
or (11) should have a common real solution. 

2 . SPECIAL CASES^ 


For n = 2 

■ 

0)^ - Kg = 0, K^aJ= 0 
file reciuired conditions are 

= 0, Kg > 0 


For n — 3- 


up « Kg 0):= 0, - Kg = 0 


fhe reauired conditions 

"K 


Eg >0, 


, ^-1 \ 
4>o. h = 4 


F or n 


4 

_ Kga* + £4 = 0, K3 0 


( 12 ) 

(13) 

(14) 

(15) 

( 16 ) 



J)iie ie(j_uir9d candit/ions arG 


,7^ > 0, k: 


- Kj^KgKj + = 0 


(17) 


Por n = 5 


5 3 A 

0) - KgO) + 0)= 0, K^o, - K £0^+ 


.gu, T 


(18) 


The required conditions are 


K^kg-Zg 


■nr^Tr” ^ o 


(19a) 


(K^K^-Kg) = (K3K^-KgK^)(E^K2- K^) 


(19tl) 


for n w 6 


6 T/- 4 __ 2 _ * tt- 5 „ 3 

(0 —Kg to -hK^to •-•Kg— 0 j — Kg 0) Kg i 


(20) 


Thece p:ive 


(ZgEg+Zg-K^K^) + ZgZg^K^Zg = 0 


+ Kg = 0 


:g(KgE:5-K^Kg)-KgiK2Kg+Zg-.Z^Ji4; 


K^(KgKg+Kg-K.,K^y-K| 



The required conditions are 


13r 


k,k/3 + k,k;_k|k,.ic§ ' “ 

and 

(K^KgKg-K^Kg-KgKg)^^ ^ 

( -iKj+K^KgKg+K ^Kg~K^ ) 

For n = 7 


(21a) 


(gib) 


7 5 3 

(0 - ICg 0 ) + K^co - Kg 0 ) = 0 (22a) 

K^to^- KgW^. KgO)^ - K^ =0 (22h) 

These give 

(JC„+K^Kg)/+ (Kg4C^K^)<fl^^(K^KgIC^) = 0 (23a) 

and 

- / [t^KgK^+K^KgK^-KjK^ 3 

+ (KgK^Kg-KgKg) = 0 (23h) 

The required conditions are 

KgK^K5~K |KgK^ +K ^K^KyKgK,^K^K2Ky+KgKgKyfK^KgKgKg->KgK6 ^ ^ 

KsV^3Vi^4+KlKgW^l^E^rWr^-^l4K5 

2 2 

+KgKgKg»“K|jIi-^+2K^K^Kg 

(24a) 


and 



136 


= (IC^KgKgKg-K^KgKg-Kll^+as^KgK^K^^^^^ 

(-KjjK g[C i+^-3^<^y-K3K4+K^K;gKgK^+K^E;gKg^^^ 

ac|+2IC^K^Kg-.IC^K|^ g+KgKgKg-K^K^) ( 24 :Td) 


Thus we find that in. every case vre get two ccaiditions. 

One of vrhich is an inequality and other is an equation* 

The second condition involves (®i~l) 

of the parameters are kept fixed, this will determine the 
value of the mth parameter. If the value of this parameter 
is changed slightly from this value, the nature of the 
equilibrium viill change, although the first condition will 
continue to be satisfied. 


S *3 TIMBl>. ] ]ELAY PIlEMTCR>;PJlSY MQhESL S 

The prey-predator model where interaction among them 
depend upon the current values is given by 


i 


^ = f(K,P) 

g = g(K.J) 


(25) 


whei-e H,S are the nimhers of prey predators respectively. 
Srowth rate of predators wlU depend not only on what they 
eat to-day hut also on the past history for the growth of 
predators* then H in the second equation of ( 26 ) can he 

replaced by 
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Q(i:) = J N(T) G(t-!I))d3} 


( 26 ) 


v/here the kernal function &(t-!E) gives the effect of eating 
one pi’oy by predator at time T after the time interval (t-l) . 
'fids gives an integrodifferential equation. Time-delay 
predator' -prey models have been discussed by Wangersky and 
anmingham and MacDonald [1976,1973 in particular . 

To solve this integrodifferential equation MacDonald 0-973 


conr.ldored the special case 


. r r-1 -az trrr 

Ur ( z) = 9. Z 6 . ) 


( 27 ) 


Q,{t) = / 




( 28 ) 


Dif-forcntiating equation ( 28 ) with respect to time 
at *“ 


= a(Q2''Q]^) 

similarly we can write 


( 29 ) 


^ = a(Q. 1-* Q^) 


( 30 ) 


g (t) = S U(®) a e 


-a(t-T) 


Now 


( 31 ) 
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on dii’feruntiating v/ith respect to time 

= / N(T) a^e~^^^“^W'+al)r(t) 

^CO 

= a(l;l-(1^) (32) 

In tliii) way one may write an ordinary differential equation 
of r'^’^ order for Q, so that the equations (25), are replaced 
"by system of differential equations given below 


m 

dt 


= f(N,P) 


^ = g(Qq»^) 

do. 

ST' ® 


d( 


j 




writing procisely the equation (33) 




e]sr(i - |) - dEP 


g = »YP + PPQ. 


dt 

dQi 

dt 

dQ, 




ar = 

Mol-cing the substitution 

H = Kn, P = K.Pj Q = Kq > T - S 


(33) 


(34) 


(35) 
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v/e ,1^0 1 


^ — bnp 

p = -cp + apcL^ 

4i= ii) 

‘ir= 

v/lipru 


b = 



J 0 




f = 


a 

e 


HI:© ma.tr ix A for this system is gi^en "by 


~ir» 

~bn* 

0 

0 

0 

« « * « 

0 

0 “ 

0 

0 

dp* 

0 

0 

* ♦ « • 

0 

0 

0 

0 

-f 

f 

0 

• « • • 

0 

0 

0 

0 

0 

-f 

f 

« * « • 

0 

0 

0 

« • 

0 

• • 

0 

• • 

0 

« • 

-f 

« ♦ • t 

0 

0 

0 

0 

0 

0 

0 

* « ■ • 

-f 

f 

^^f 

0 

0 

0 

0 

• • • « 

0 

-f 


(36) 


(37) 


(38) 


I’he eciuatim for eigenvalues on sinjjlificat ion, gives 

j^r+2 ^ ^r+l ^ ^ ^r 

+ n* ^Ogf^) + ... 

+ (f^+n* + Xn*f^ + af® = 0 (39) 


n*=cL* a 






fr‘» “ == I 


(40) 


whoro 
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v/e get 

li = n(l-n) - bnp 
p = _cp + dpq.^ 

4i= <ip 

q^= f(3a-q^) 

where 



The matrix A for this system is given by 


~n^ 

-bn* 

0 

0 

0 

« • • * 0 

0 ~ 

0 

0 

dp* 

0 

0 

% * m » 0 

0 

0 

0 

-f 

f 

0 

• « • • 0 

0 

0 

0 

0 

-f 

f 

.... 0 

0 

0 

0 

0 

0 

-f 

.... 0 

0 

0 

0 

0 

0 

0 

.... -f 

f 


0 

0 

0 

0 

.... 0 

-f 


(36) 


(37) 


(38) 


Tbfi eq,uation for eigenvalues cn simplification, gives 
^r+2 ^r+1 ^r^^^ ^r 

+ n* ^Ggf^) + ... 

+ (f^+n* + xn*:? + af^ = 0 (39) 

JL 


where 


n*=q^ 

1 


a-» = 
^ 2 


.= = 


d’ 


r, 

J* _ 


d-c 

W 


a = 


(d^c) (40) 
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5.4 

SPEGIAl OASES 


(i) 

11 



(f + n*) + xn*f 4- af = 0 

(41) 

(ii) 

r = 2 



x"^ + x^(2f+n*) + x^(f^+2n*f )+ xn*f^+af^= 0 

(42) 

(iii) 

r = 3 

X^ + x^(3f+n*) + X^(3f^+3n*f) + X^(f^+3n*f^) 



+ xn*f^+af^ = 0 

(43) 

(iv) 

r = 4 



X^ + x^(4f+n*) + x^(6f^+4n*f) + x^(4f^+6n*f^ 

) 


+ X^(f\4n*f^) + Xn*f\ af^ = 0 

(44) 

Por discussing the existence of points of bifurcation, 

we 

apply the conditions (15), (17), (19) and (2l). 

Por r = 1 the required condition is 



d^-dc-c 

1 _ — 

(45) 

Thus 

for a bifurcation point to exist, we should have 



d^ > (d+l)o 

(46) 

!■ or 

r=2, the required condition is 



-2n*f^+f^(4a~4n*^)+f(4a%.-2n*^) + an*^ = 0 

(47) . 


This equation does have a positive real root and as such 
in this case, a hifur cation point does exist. 
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i'or r = 3, the required conditicaas are 

+ n*) - ^ Q 

(3f+n^)(3f2+3n*f) - (f3+3n^^f2) 

and C^-t -vO 

~8ii^f\f^(E4a-24n^^^) +/(-.24n*^+ 45 an*) 
+f (-.8n*^4-34an*^+a^) + 9ocn*^ = 0 


(48) 


(49) 


Equaticai (49) has at least one real positive root* If this 
satisfies oqn. (48), a point of hifurcation exists. 


I'or r=4, the required conditions are 

( 4f +n* ) ( 6f ^+4n*f j (n*f ^)~ ( 4f+n''®' )^af ( 4f ^+6n*f ^ )n*f ^ 

(4f+n*) ( 6f ^+4n*f ) (4f^+6n*f^ )+(4f+n*)n*f^-(4f4-n*)^(f^+4n*f^ 

_(4f^+6n*f^)^X50) 

and ■ 

-.4ai*^f^^-.f^°(n*\480n*^) + f^(32n*^+ 576an*-824n*^) 
+f®(ien*^-528n*\704n*^a-256a^-.1280a+320n*) 

+f ’^( 8 16n*^- 152n*^+ 128n*^~64n*-5 184ocn*+3 72an*^ 

+256a+2n*^-256n*a^)+f *^(65n*%632n*^ 

-.74‘24an*\l232n*^+576cm*+90an*'^-.l6n*^-.96n*^a^ ) ) 

j_ ^ 3 5 2 

+f °(208n* +2156h* +8n**^236dai* +I024an* 

-.2752cai*\8n*^«-16n^V)+f^(24n^+2527n*4 

-384n *\23o4an*^~ 1344ocn*'^«n*\^+'?6hn ^‘^ ) 

+f^( 832n*^-96n*4-1024a*^- 144001*^+ 1344an*^+6an*'^ } 

+f^(96n*^-768n^^+ 304c£n«4 -8n^) + f (24cm*^-192n*^) 

(51) 


16n*^ = 0. 
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Equation (51) will haw no real positive root or it will 
iiaA?e even nuEitier of real positive roots. If these saitisfy 

equation. (5o),; on even number of a points of bifurcation will exist. 


5.5 BIEUROATIOW TIIE50RY ECR THE GROWTH OP PaPUIATIOH 


Pollov/ing May (1973), we tafee 

ij 

S = I®~bH^~GN / H(t» )G(-b-t» )<it' (52) 

^OO 

the Icernal function G(t-t*) gives the influence of the 
population at time t’ on the populatim at time t. 


Let 

and 


-.a(t~t’ ), , sP~l p 
G(t-t*) = — -a- 


Qp(t) 


t -a(t-t’ ) / . . , R-1 

= / N(t«) 

r&) 


at' 




then we have 

~ = KH(t) - bH^(t) - C[H(t) Qp(t) 
dQ 

^ ^ a(Q^(t)-Q2(t)) 

^ = a(H(t)-Q^(t)) (53) 
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S' or equili'brium 


IT = ^ ^ Ql = Qg 


K 


• • • 


The matrix A for the system Is ^siven by 

-hi -.gTt 0 

0 -a a 

0 u -a 


0 

a 


0 

» ft 

ft ft 

0 

0 


0 

0 

a 


0 ... 0 

0 ... 0 

0 » * « 0 


0 

0 


0 

0 


0 

0 


-a 

0 


0 

0 

0 


. . 

a 


(54) 


The equaticn for eigenvalues, on simplification, gives 


(^'c^a + hi) + x^“’^(^Oga^+^O^ablT) 

+ X^ ^ Gga^bN )+ . . . 

+ bla^+ a^Ct *= 0 (55) 

5 -5 SP ECIAL OASES 

(i) P = 1 

X^+ xCa+bl) + a^(b+C) = 0 (56) 

(ii) P = 2 

X^4. X^(2aH-bi) + + a^(M+ct5 = 0 (57) 

(iii) p = 3 

x\ X^(3a+M) + X^C3a^+3abN) 4- X (a^Sa^#) 

4 a^CblT 4 ®) 0 


(58) 
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(iv) p = 4 

X^+x'^(4a+bN)+x^(6a^+4aM) + X^(4a^+6a^bir) 

+ X {a^+ 4 bia^) + a%(b+C) = 0 (59) 

lor discussing the points of bifurcation, we apply conditions 
(13), (15), (17) and (I 9 ) 

lor p = 1 

a+bK = 0, a|[(b+c) >0 ’ ( 60 ) 

there is no real positive value of a so there ?7ill be no point 
of bifurcation, 

lor P = 2 

2 a^+a( 4 bif-aN) + 2b^i^ =0 ( 61 ) 

Tliis equation either has no positive real root or gives even 
number of real positive roots. If it gives even number of 
positive real roots, there exist two points of bifurcation. 

lor p = 3 

>0 ( 62 ) 

3c!H-bIT 

and 

(a^+SaSi)^- (3a+bl) (3a^+3abN) (a^Sa^blT) 

H- a^(bl+(ir)(3a+bl)^ = 0 (63) 

~8a^+a^(90~24b)+aK^(6bG-24b^) + i^(b^C-8b^) = 0 

(64) 

If C/b > 8, point of bifurcation exis'isif that value of a also 
satisfies eqn, (62), 
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If C/Td < 8, then either there will escist no positive real 
root or there exist even number of positive real roots. 

If roots exist and satisfy eqn. (62) there will exist a point 
of bifurcation. 

lor p = 4 the required conditions are 

(a\4b&a^)(4a+bl) ^ a%(b+Q) ^ ^ 

(4a+bN) (6a2+4abH)-(4a3+6a2b^) 

and 

64a\ifa^ ( 256b- 1120 )+a^^ ( 384b 144bC-0^ ) 

+aH^(256b^-8aib^C)+ir^(64b^-16b^O) = 0 

If 4 < 0/b there exisi:3a real positive root. If this 
satisfies eqn, (65) there exist a point of bifurcation. 

If 4 > 0/b then either there will exist no positive real 
root or there exist even number of positive real roots. 

If roots exist and satisfy eqn, (65) there vriH exist a point 
of bifurcation. 


(65) 


( 66 ) 



CHAPTER 6 


TRMSITIOI EROM DETERMBTISTIG TO STOCHASTIC MODELS 

6.1 DETERIvIIHISTIC M’D STOCHASTIC MODELS 

All the models considered so far isii the thesis, 
whether in terms of differential equations or difference 
eq.uat ions or integro- differential equations or delay- 
differential equations, have heen deterministic i.e. our 
object has been to deterniine the population sizes of the 
various species as definite and uniquely defined functions 
of time. In these deterministic models for population 
growth, Y'je seek to determine the population size H(t) as a 
function of time. If the model is assumed correct, there 
is no uncertainity about H(t). 

In practice however chance plays an important role and 
very often v/e prefer to seek to determine p(n,t) ^/hich is 
the probability of the population size being n at a given time 
t. Thus at a given time t, the population size does not nave 
a fixed value but may take values 0 , 1,2,3, .. .,n. . . with 
probabilities p(0,t), p(l,t), p(2,t), p( 35 t), . . . ,p(n,t) ... 
so that at any given time we seek to determine the probability 
distribution of the population size. We, of course get a 
family of probability distributions one probability distribution 
for each t. Such models for determining probability distributions 


3 X 0 CcillBd sijOC ho-S’b ic modBls 
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In the present chapter, we develop a general method 
of f ormula 0 ing the stochastic models, corresponding to a 
deterministic model given in terms of a system of differential 
equations. Ihis method v/as first initiated by Kapur pL97i'j , 

The basic idea is to regard the positive terms in the 
expressions for the derivatives of population size as due to 
‘births* and the negative terms there as due to ’deaths*, so 
that every system corresponds to a non-linear birth and death 
process. This approach gives in general an infinite system of 
different ial-d iff erence equations for determining the infinite 
set of probabilities p(n,t). In general the differential 
equation determining p(n,t) involves p(n+l,t) and so if tadce 
a finite subset of equations, it cannot determine all the 
probabilities. As such v/e have to consider the infinite system 
of differential difference equation as a whole. 

To solve the system v/e therefore take recourse to the 
use of probability generating fijnction defined by 

CO 

<{.(s,t) = Z p(n,t) s^ 
n=0 

We determine a partial differential equation for solving for 
(|,(s,t). If we can solve this equation, we can find <^(s,t) 
and then by expanding it in power of s, we can find p(n, c) . 

This method has been successfully applied to linear birth- 
death— immigration processes. However when the processes is 
non-linear or when emigration is involved, some difficulties 
arise in solving the resulting partial differential equation. 
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Por linear Id irtii-death- immigration, emigration process, 
a solution was given by G-etz Qi975] . However he omitted a 
term containing the probability of extinction and thus all his 
calculations were wrong » Later Kapur pointed out this mistalce 
which was accepted by G-etz follovying cases wore 

discussed by Kapur. 

(i) Birth-Leath- emigration process : Ihe complete solution 
¥/as obtained in terms of incomplete beta function 

Kap ur j”l9 79 Cl , [19 79 . 

(ii) Birth-Leatli-emigrat ion process with birth rate less 
than death rate. Ihe solution was obtain by solving 
an Abel's integral equation (Kapui‘ and Saleem [197^). 

(iii) Birth-Deatluiimnigration-emigration process : Ihe 
complete solution for the steady case was obtain oy 
Kapur [;i978a, 1978b, 1979c, 1979d3 and (Kapur and 
Kapur [19 78] ). 

Ihe processes which arise in general a,re non-line an 
and it is not possible to solve the partial differential equation 
for the probability generating function for the following reasons, 
(a) Ihe partial differential equation contains in addition 

to the unknovn function <jl(s,t) some terms involving 
unknown probabilities. 

'Ihe partial differential equations are highly non-linear. 
v7e have the initial constrains in the form 

(|> (s,0) = S ° 


(^) 

(c) 
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where is the initial number of person in the system, but 
no boundary conditions are given so that for the steady case 
we have to solve partial differential equation without any 
initial and boundary conditions. 

We have investigated the partial differential equations 
which arise in a number of problems. We have not attempted to 
solve these. However we have obtained with their help 
differential equations for moments for the probability distribu- 
tions, Unfortunately aigain here the system becomes indeterminate 
in the sense that the number of unlmown moraents to be found out 
are always more than the number of equations. We can hov/ever 
integrate these in some cases to obtain some interesting 
relations betv^een the moments. 

6 .2 I'ORI'iAflOH OF IHS STOGIhiSlIG MODS! GORlSSPONUn^iG TO A CtIYM 
DHl'EfllBllSflG IIOUEL 


Gonsider a general deterministic model 



n 

2 

0=1 


a. 


10 





n 

2 

0=1 


b . . X . + c . 

10 0 1 


( i — 1,2, 


n) 

( 1 ) 


Here all the coefficients are positive, let p(X^,X2, , , .,X^, t) 
be the probability that at time t, there are X^ individuals 
of the first kind, Xo individuals of the second kind and so on. 
Also we assume 


P (birth of an individual of the ith species in time At) 


= ( X. + 0^) At + o|At;^ 
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P (deatii of ail individual o-i tlie ith species in tine At) 


= ( S X. Ij ^t + 

1 J 


P (more than one Ghenge in time At = 0(At)) 
so that pro'babilit 3 r of no change in time At 


n n 


n n 


= 1 - ( S 1 a. .X, + E c J At -( 2 Z t..X.Xj At 
i=l 3=1 J j=l * i=l 3=1 - 3 


+ or(At)j 


We then get. 


n n 


p(X^,X 2 ,..M-C^?t + At) = p(XpX2,...,2^,t) Z 1«( Z 2^ a. . X. 

i=l 3“-^ 


+ Z *^ •? ) At 
3^=1 ^ 


n n 

( Z S b X I 
i=l 3=1 ^ J 


,) ^' 1=3 


+ .1 p(Xi.Xb ••• ^>■‘5 

1-1 

r. A ^li + ^ii(Aik) + nH 


3=1 ^ 

D 5 ^i 


+ Z p(X^,ZCg, . . ^x+1 ••• 

D -1 


C .2 X. (X.+^) -At + 0lAt,3 

3=^1 
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Pi" oceGding co "bhe limit; as “*■ Oy W6 gsd 


r ^ ^ 


i=l 0=1 


ID D 


+ I 0. 


n n 

- p(X.,Xp,...yX ,t) s S b.. X. X. 
^ ^ ^ i=l 3=1 ^ ^ 




C jfl ^iD ^3 °i ^ ^il 


4- H ^ p Xg 9 • • • ^X^-4* 1 j 9 ^ 


J 4- 


n n 


= -p(X^,Di2,...,Xri’^^ C .2 .2 


i?=l 3=1 ^ 


2 °i3 


n n 


- p(X2_,Xg, ^2^ 


wL* 

+ I pCX^yXgy ...,Vl’^i+1» 

i=l 

c A ^ij - V + “iD 
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n n 


3=1 


+ i>ii (Xi*i) :] 


(4) 


These eq.uations will hold for values of S^,2gj . , . ,2^ frou 1 to 
with slight modifications, these will also hold for the cases 
when some X’s are zero, 

How we define the generating function 

00 00 C50 

$ ( s -1 ) Sp j • • • » ®Yi j t ) = 2 2 2 p (X-| , Xp , • , • j}i,^ j t ) • 

X^=0 Xg^O x^=o 


^2 \ 

• "2 ••• ^n 


(5) 


X. Xg X 


Multiplyuig (4) by . Sg^ ... and summing up for all 

X’s, we get 


|t = - “ S a . s.(l-sO l|--+ 2; .S, 

3=1 is=l 3 3. 00 ^ 13 J J- 


8<}>, 


n n 


6 *t’ 


3 3=1 j=l 


n 


fel 


+ _s. ty. Oj, (L.3,) 


n 


( 6 ) 


i=l 


+ some terms independent of (|) , 

If there are X°, 2^, ...,X^ individuals at time t = 0, we get 
■the initial condition 


0 

1 


4> (s^,Sg, . ^ . ,s^,0) = Sg 




C?) 



153 


6.3 GOEIPETITIOJNi FOR FCOE ^\K) SPECIIS 


H = - bXI 

dY 

= PY - clY 


( 8 ) 


In the absence of other species, each species grov/s exponentially 
and both species complete for the same food. Hence all the 
coefficients are positivre. 

Lot p(Y,Y,t) be the probability that at time t, there 
are X individuals for the first Rind and Y individuals for the 
second kind. Also we assume 
Birth of the first species at time t = aX 
Death of the first species at time t = bXY 
Birth of the second species a;t time t = PY 
Death of the second species at time t = cU 

Probability of more tnan one change in time At = 0(At) we then 

get, 

p(X,Y,t+At) = p(X-l,Y,t) [;a(X~l) At + c[(Atj|] 

+ p(X+l,Y,t) ];b(X+l)Y At + 0(At)3 
+ p(X,Y-l,t) j5(Y-l) At + O^At)j]] + p(X,Y+l , t) 

[;;cX(Y+l) At + C^(At>3 + p(X,Y,t) Cl^CaX+bXY 
+ PY + 0X1)3 At + ci(At)-': 


( 9 ) 
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Proceedmg to the limit At - o, we get 

= p ( a— I j Yy t ) (l— l)a + p(Z+l5Yjt) (X+l) Yb 


+ p(lyY-l,t)(Y-.l) P + p(X,Y+l,t) X(Y+l)c 


- p(X,Y,t) [^aX + b XY + PY + cXY”| (10) 

X'^. 

Similarly v/hen the population of first species "’is zero, vve get 

p(0,Y,t) = p(l5Y,t) Y b + p(0,Y-l,t) (Y-1) P - P(0,Y,t) IP 

( 11 ) 

\;hen the population of the second species is zero 


—[j p(X,0,t) = p(X-l,0,^i:)(X-l)a + p(X,l,t) Xc - p(X,0,-t) Xa 

— (12) 

and v/hen populations of both species are zero at time t mq get 


p(0,0,t) = 0 (13) 


Summing (lO), (H)? (12) and (13), we get 


oo oo oo oo 

|-p(X,Y,t)= 2 Z p(X~l,Y,t) a/(X-.l) + Z Z p(X+l,Y, t)(X+l)Yl 

X=1 Y=1 X;=l Y=1 

/ 

oo oo ' ^ oo oo 

+ Z Z p(X,Y-l,t)i(Y-l)P+ Z Z p(X,Y+l,t)X(Y+l)c 
3&=1 Y=1 1^1 Y=1 

/ 

// 

OO oo '• 

- Z Z p(X,Y,t^) (aX+PY+bXr + cXY) 

X=1 Y=:l / ■ 

+ Z p(l,Y,t) ;bY + Z p(0,Y-l,t)(Y-l)P 
Y=1 Y=1 
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~ 2 p(0,I,t) IP + Z p(X-l,0,t) (X„l) 

X= 1 -1 




a 


+ 2 p(X,l,t) Xc - Z p(X,0,t) Xa 


X=1 


i 


(14) 


Now we define "the genera.ting function 


oo oo 


<J>( 2 ;,s,t) = Z 2 p(X,Y,t) 
iL?=0 Y=30 


(15) 




Vi 


Multiplying equation (1^) thi-oughout by z^ s^Aand summing for 
all X and I, we get with the help of (15) 


|4 = e(b-i)s h* 0. (Us) + 


(^1) H 


(16) 


6.4 LAUNCHES TIR MODEL I'OR A BATTLE 


d?b’ 


= P - al 


(17) 


dl 

re 


= q - bX 


where X and I are the numbers of effective attacker and defender 
troop, P and q are the rates at v/hich thej’ are being rein.forced 
and a and b are attacker and defender casual ities per opposing 
combatant. 


Resulting partial differential equaxion by taking 
generating function 


“ Y r 

<t>(z,s,t) =22 p(X,I,t) 2"" s'" 

X =0 1=0 


( 18 ) 
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is given by 


- 

r ^ a . 

• 3«S 

J. 

9s j 



as 


a 


2 Y p(0,Yj-&) 
1=1 


.Y 


(1-1)' 


+ b 


2 X p(X,0,t) (1 - 1 

^1 ^ J 


1^1 

+ q.4) (s-l) +'^(|) (z-1) 

6.5 MODEL I'OR lEYERLAilOIAL OBIADE ^ 
dX 


dt 


= -aZ + bZY 


dY 

§ = -PY H- err 


(19) 


( 20 ) 


Here trade is beneficial to both nations and isolation is 


harmful. 


Resulting partial differential equation by talcing 
generating function 

~ 7 Y 

({((ZjSjt) z= E 2 p(X,Y,t) z"'*' s 

X=0 Y=0 


( 21 ) 


is given by 


5 <f' 


as 


P(l-s) 






1 . 
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6.6 PREDATOR-PREY MODEL 

c\Y 

^ = ax ^ bXY 
dY 

§ = oXI - PY 


(23) 


Here S, Y are the two populations at time t. In the absence 
of Y species the species'! grows exponentially while in the 
absence of X. species, the species Y decreases exponentially. 

The mutual contacts are useful for Y and harmful for X species. 

By talcing generating function 

oo oo 

(|>(z,s,t) =21 p(X,Y,t) '/• s^ (24) 

2=0 Y=0 

the resulting partial differential equation is 


6<() _ 6^ 
rb’ ~ as 




+ H Caz (z-l)2l 


+ hh s(s-l) + bs (l-2)n 

6.7 MUiBlESHIP OP POLITICAL PARTIES 
= a2 •“ bX + c 

H = -P2 + q.Y + r 


(25) 


(26) 


Here the membership of a political party has an encouraging 
effect on its own strength and the membership of the opposite 
political party has an inhibiting effect. 
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Resulting partial differential equation by taicing 
generating function 


oo oo 

<|)( 2 ,s,t) = i: I p(X,yyt) s^ (27) 

X=0 Y=0 

is given by 


If = ll-lMCs-l) + te(|- 1)3 

+ <|)[Z^(s-l) + c(z-l)21 + ^ 

oo 

+ P 2 Z p(x,0,t) (1 _ 

Z=1 


+ || Caz (z~l) + 

OO 

2 I p(0,Y,t) s^ 
Y=1 


P_z 

s*” 

(1 


- Pz 


a 



(28) 


6 .8 ARMS RACE BE1\VEEM NATIONS 


= — cllL -h "bX -f C 
S = H-<lY+r 


(89) 


The rate at which a nation build up its armed-strength 
is proportional to the armed strength of the enemy and also 
to its own, the rate is positively correlated with en(Sii 3 ’' 
strength and negatively cor relented ’.vith its ov?n. 

The resulting partial differential equation by tailing 
generat.ing function 

<|>(z,s,t) = 2 2 p(Z,Y,t) z^ s'" (30) 

X=0 Y=0 


is given by 
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S<j> 

aT = 


M 

ds" 


Ila(l-s) + bs(z-.l)^ + |1 [;pz (^_i) + aiUz)2 


+ <})[]c(z-l) + r(s-l)[] 


(31) 


6.9 MOMMTS 01 PR03i\BILHT OIO'IRIBUIION 


lo deduce the moments, v/e use the relations 


m) = 11 Ki.) = H- E(xr) = O 


5 zes 


(32) 


_ LU j. Li v{'^\ , s 


= §4 + |i .(|i)t Var(X) 4-4 ^ H -ilif 


(33) 


where all the partial derivatiTes are evaluated at a = 1, z = 1. 
lor any model, to obtain the moments we differentiate the 
partial differential equation of the probability generating 
function for the model. 


(i) Once Y/ith respect to s, once with respect to z and 
then put s = 1, z = 1. 

(ii) fwice with respect to s, twice \ 7 ith respect to z, once 
vTith respect to a and then with respect to z and then 
put s = z = 1 and so on. 

Using this method we get the following relations 
(a) AtaiS S.AGE lA'IIOHS 

^ B(X) = b E(T) - a E(X) + c (34) 



160 


dt + PS(iv) + r ^ 35 ^ 

^ E(X^) = aE(X) + bE(I) + 2c E(l) + 2b E(ir) 

- 2 a E(X^) + c ( 36 ) 

E(Y ) = q.E(Y) + 2rS(Y) + PE(X) + 2PE(XlO 

-2q.E(Y^)+r ( 57 ) 

E(XY) = bE(Y^) + PE(X^) - (a+q) E(XY) 

+ cB(Y) + rB(X) (38) 

(b) HvURlLVIIOML 'Iti^iES 

^ E(X) = -a E(X) + b E(XY) (39) 

E(Y) = _P E(Y) + c E(XY) (40) 

^ E(X^) = a E(X) + ,^b E(XY) 2 a B(:^) +^bE(X^) ( 4 l) 

|:g- ECy"") = PE(Y) + 2 rc E(rf) - 2 P E(Y^) +^c E(XY^) (42) 

E(XY) = -(P+a) E(XY) + c E(X^Y) + bE(XY^) (43) 

(c) PRIIDAI'OR - PRBI MODEL 

^ E(X) = aE(X) - bE(XY) (44) 

l=E = -PE(Y) + cE(XY) (45) 

S(X^) = aE(X) + bE(XY) + 2aE(X^) - 2bE(X^Y) (46) 
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E(I^) = PE(Y) + cE(XI) - 2PE(Y^) + 2cE(Y^Z) 

E(XY) = (a^) E(XJ) + cE(j^Y) - bS(ZY^) 

(d) Gasp EDITION MODEL 

E(X) = aE(X) ~ •bE(XI) 

E(Y) = PE(Y) - CE(XY) 

E(X^) = aE(X) + 2aE(X^) + TdE(IY) ~ 2bE(X^Y) 
^ E(Y^) = PE(Y) + HPE(Y^) + cE(XI) - 2cE(:CI^) 

E(XY) = (P+a) E(XI) - bE(XY^) - cE(X^Y) 

(e) LADNOHESTEt MODEL 

E(X) = -aE(Y) -1- a 2 I p(o,Y,t) + P 

Y=1 

oo 

^ E(Y) = -bE(X) + b 2 X p(X,0,t) + q 
^ X=1 

OO 

E(r) = aE(Y) - a 2 Y p(0,Y,t) 
a-c Y=i 

+ 2PE(X) - 2aE(XY) + P 

CO 

|r E(Y^) = bE(X) - b 2 X p(X,0,t) 

:&=1 

+2q E(Y) - 2bE(XY) + q 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 



162 


3 ^ E(Xr) = _aJ3(Y®) _ tBtX®) + (iDd) 

<50 ^ 

p oo 

+ PB{X) + a D 1“= P(0,l,t) + b Z X® P(2,0,t) 

1=1 jr_2 

(58) 

(f) KMBmS OP POLIPIG.il PilPflES 

B(X) = aE(S) - bE(I) + b 2 Y p(0,I,t) + c (59) 

Y=1 ^ 

It + P 2 X p(X,0,t) + r (60) 

X=1 

d 2 

S(X ) = (2c+a) E(X) + bE(Y) - 2bE(XY) 

2 “ 

+ 2aE(X ) - b 2 Y p(0,I,t) + c (6l) 

1=1 . 

^ S(Y^) = (q+ 2r) E(Y) + PE(X) + 2qB(Y^) 

OO 

“ 2PE(XI) - P 2 X p(x,0,x) + r (62) 

2>=1 

E(ia:) = (cn^ E(]CY) - PE(Z^) ~ bE(Y^) 

+ + rE(X) + cECY) + b 2 Y^ p( 0 ,Y,t) 

Y=1 

^ 2 

+ V E X p(X,0,t) (65) 

X=1 

6.l0 RELAl’IOlI BET\/EM MOBIEtiPS 
(i) ARMS RACE BETiffiM 'HW RATIONS 

We choose i and m positive constijits to satisfy 


relation (64) 
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injp^ - aji ab - mq. 

then we get from equ8.tion (29) 

^ C^iE(X) + mE(Y)^J = k [; ilE(X) + mS(T)3 



+ £c + mr 

(65) 

£ = (q_a) + f(q^a)^ 

+ 4bP, m = 2b 

(66) 

■f(q-a)"^ + 4biP - 

k-i-a) 

(67) 

It - 2 



On integrating equation (65), we get 

JlE(X) + mE(T) = ilX^ + mY^;] (e^^ - 1) (68) 

v;here X^ and Y^ are the initial values of X and Y. 

Equation (66) sliox^ that 1 !E(X) + mS(Y) will increase or 
decrease as empt otic ally with time according as 

hP I aq (69) 

We may note that for the deterministic model 

, . £,c +- mr , , , 

jlX + mY = (illg + biYq) e " + (e - 1) (*70) 

% 

so that ax + mY for the deterministic model and ilE(X) + ml(Y) 
for the stochastic model are identical function of time. Ihis 
is obvious otherwise since X, Y and E(X), E(Y) satisfy the 
same differential equations. In this cose from (34) and (35), 
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we can find S(X) and E(Y) as function of time and then can 
oi^tani E(X^), E(Y^), E(XY) from (36), (37) and (38). 

(J-i) POR miEEUUTIOHAL 'liRALE 

_ ^ vAv-y ' ^ ' 

Prom Sq.uation (.SO)', we get 

= ^PE(Y) ~ acE(X) (71) 

This corresponds to tiie equation 

(cX - bY) = bPY - acx (72) 

for the deterministic model. 

If hov'/ever a = P, v/e get 

[lcE(X) - bE(Y);] = -ai:cE(X) - bE(Y)3 (73) 

so that 

cS(X) - bE(Y) = (cXq - bY^) (74) 

and cE(X) - bE(Y) decreases exponentially v/ith time and 

It C *- ^ 

t -* oo 

f i ii ) EOS PKEDATOR..PREI MOEBE 

Prom equation (■S3'), we get 

[;cE(X) + bE(Y)3 = acE(X) - bEE(Y) (76) 

Phis corresponds to the equation 

Cox + m 3 = aoz - 1*Y , C77) 


for the deterministic model. 
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(iv) FOR OOI,IPE1'ITI01\T MOREL 

From eOLuation (S'), we get 

C;cS(X) - bE(Y)3] = acE(X) - bPE(Y) (78) 

Lliis corresponds to 

QcX - bY^ = acX bPY (79) 

for the deterministic model. 

If a = P , ?;e get 

^ - bE(Y);] = a[;cE(X) - bE(Y)3 (8o) 

so that 

gE(X) - bl(Y) = (cX^ - bY^) (81) 

Thus cE(X) - bS(Y) increases exponentially with time unless 

"^0 = ^""o 

(a) If cXq = bY^ then cE(X) = bE(Y) for all t (82) 

(b) If cXq > bY^ then cE(X) > bE(Y) for all t (83) 

and 

cE(X) - bE(Y) - ~ and t ~ (84) 

In this case for the deterministic model X - Y - 0. 

(c) If cXq < bY^ then oE(X) < bE(Y) for all t and 

cE(X) - bE(Y) and t - “ (85) 

In this case for the deterministic model X - 0, Y - 
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